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1.0  RESEARCH  OBJECTIVE 


Impinging  jet  injectors  are  a  common  type  of  injector  used  in  liquid  rocket  engines  and 
are  typically  used  in  engines  where  both  propellants  are  injected  as  a  liquid,  e.g.,  engines  using 
LOX/hydrocarbon  and  storable  propellant  combinations.  Because  they  are  much  simpler  to 
fabricate  than  coaxial  injectors,  they  are  also  an  attractive  candidate  for  use  in  LOX/hydrogen 
engines.  There  are  also  current  interests  in  engines  where  both  propellants  are  injected  as  gas 
and  where  three  propellants  are  used;  here,  the  impinging  jet  injector  system  would  also  be 
favored  over  coaxial  jets  due  to  superior  mixing  characteristics  and  potentially  simpler 
manifolding  arrangement 

The  current  state  of  impinging  jet  injector  design  analysis,  however,  significantly  lags 
behind  that  for  coaxial  elements.  A  substantial  problem  exists  with  the  current  design  analysis 
methodology  for  impinging  jet  injectors  in  that  to  increase  the  margin  from  unstable  combustor 
operation,  the  combustor  designer  is  faced  with  decrementing  performance  and  thermal 
compatibility  characteristics.  Also,  empirically-based  analysis  techniques  are  used  which 
depend,  of  course,  on  existing  engine  data.  Due  to  the  empirical  nature  of  stability  analysis,  and 
the  very  high  costs  of  engine  development,  innovative  designs  that  take  advantage  of  current 
technology  are  not  being  used;  instead,  current  engines  are  essentially  being  built  with  the  same 
injector  designs  that  were  developed  in  the  1960's. 

The  research  program  at  the  Propulsion 
Engineering  Research  Center  at  Penn  State  has 
been  focused  on  providing  the  requisite 
fundamental  understanding  for  the  development  of 
an  advanced,  a  priori  combustion  stability  design 
analysis  capability.  This  can  only  be  achieved 
when  the  fundamental  mechanisms  of  combustion 
stability  in  rocket  engines  using  impinging  jet 
injectors  are  understood.  A  key  element  of  the 
program's  approach  was  the  use  of  a  proven 
empirical  stability  correlation,  previously  a 
proprietary  design  tool  at  Aerojet  shown  in  Fig.  1, 
that  relates  the  highest  possible  frequency  of 
combustion  instability  that  can  be  driven  by  an 
injector  with  given  geometric  and  operational  parameters,  i.e.,  injector  diameter  and  injection 
velocity.  A  review  and  characteristic  time  analysis  of  potential  combustion  instability 
mechanisms  that  explain  the  correlation  has  been  performed,1  resulting  in  the  identification  of 
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Fig  1.  Hewitt  stability  correlation  for  deter¬ 
mining  the  highest  possible  combustion  in¬ 
stability  frequency  in  rocket  engines  using 
impinging  jet  injectors  as  a  function  of  injec¬ 
tor  stability  parameter.  Correlation  effective¬ 
ness  is  illustrated  by  a  few  examples  of 
LOX/hydrocarbon  engines. 
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the  processes  of  primary  and  secondary  atomization,  inter-propellant  mixing,  and  vaporization  as 
likely  important  controlling  steps.  Identification  of  these  processes  was  used  to  define  the 
dependent  test  variables  to  be  measured:  length  of  the  intact  liquid  sheet,  ligament  shedding 
frequency,  drop  size  and  velocity  distribution,  and  heat  release  rate.  The  empirical  correlation 
was  used  to  define  independent  test  parameters:  orifice  diameter,  injection  velocity,  and  imposed 
frequency  of  oscillation. 

The  results  of  the  review  also  indicated  that  a  clear  understanding  of  primary  atomization 
is  of  particular  importance.  Studies  of  primary  atomization  have  shown  that  it  occurs  on  a 
temporal  scale  that  is  similar  in  magnitude  to  combustion  instability,  and  that  the  temporal 
dependency  on  orifice  diameter  and  injection  velocity  is  similar  to  the  dependency  indicated  by 
the  combustion  stability  correlation  used  in  the  present  study.  Furthermore,  inter-propellant 
mixing  and  drop  vaporization  and  burning  are  controlled  by  the  distribution  of  drop  size  and 
velocity  which  is  determined  by  primary  atomization.  Thus,  the  major  effort  has  been  to 
characterize  the  injector's  geometric  and  flow  effects  on  primary  atomization  and  to  develop  a 
reliable  and  accurate  atomization  model  for  impinging  jets. 

In  summary,  the  overall  objective  of  this  study  was  to  determine  the  physical  mechanisms 
that  are  responsible  for  the  combustion  stability  characteristics  of  rocket  engines  that  use 
impinging  jet  injectors,  using  a  proven  empirical  correlation  to  link  the  fundamental  study  to 
practical  rocket  combustors.  Additionally,  the  most  complete  and  systematic  study  of  the 
atomization  of  impinging  liquid  jets  to  date  has  been  completed  under  this  contract  Effects  of 
jet  flow  condition,  orifice  diameter,  impingement  angle,  pre-impingment  length,  fabrication 
procedure,  and  jet  velocity  at  steady  and  oscillating,  and  atmospheric-  and  high-pressure 
conditions  have  been  investigated.  From  these  cold-flow  tests,  three  important  observations 
related  to  the  basis  of  the  stability  correlation  have  been  made:  (1)  as  the  stability  parameter  is 
increased,  leading  to  more  stable  operation,  the  spray  drop  size  distribution  becomes  more 
polydisperse;  (2)  as  the  stability  parameter  is  increased,  the  mean  drop  size  is  increased;  and  (3) 
the  frequency  with  which  ligaments  are  formed  and  with  which  waves  of  drops  enter  the 
combustion  zone  has  a  functional  dependence  on  the  stability  parameter  similar  to  that  of  the 
maximum  predicted  instability  frequency  as  shown  in  Fig.  1. 

The  results  of  this  study  have  contributed  to  a  fundamental  understanding  of  the 
atomization  processes  and  the  operative  instability  mechanisms  associated  with  impinging  jet 
injectors.  These  results  must  still  be  verified  under  combus ting-flow  conditions.  Also,  more 
work  needs  to  be  done  to  relate  the  observed  phenomenology  to  specific  physical  processes  to 
enable  development  of  a  priori  design  analysis  methodologies  and  implementation  of 
engineering  control  strategies  for  future  design  of  stable,  high-performing,  and  thermally 
compatible  rocket  combustors. 
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2.0  RESULTS 


The  subject  study  was  concentrated  on  defining  the  operative  mechanisms  of  combustion 
instability  in  rocket  engines  that  use  impinging  jet  injectors.  General  information  regarding  the 
combustion  process  in  rocket  engine  combustors  of  all  types  is  also  an  important  byproduct  of 
this  research.  A  review  conducted  earlier1  identified  the  combustion  processes  of  primary 
atomization,  secondary  atomization,  inter-propellant  mixing,  and  drop  heating,  vaporization,  and 
burning  as  potential  key  mechanisms  of  combustion  instability.  The  effort  involved 
cold-flow  studies  of  atomization  under  both  steady  and  forced  oscillatory  ambient  conditions. 
There  has  also  been  a  parallel  effort  in  developing  an  accurate  atomization  model.  More  details 
of  the  work  can  be  found  in  References  1-3,  and  in  Appendix  1.  Future  studies  under  high- 
pressure,  combusting,  and  oscillatory  conditions  must  still  be  undertaken  to  verify  the  cold-flow 
results  and  to  measure  the  combustion  response  of  impinging  liquid  jet  spray  combustion  in 
terms  of  the  phase  relationship  between  the  oscillating  pressure  field  and  combustion. 

The  development  of  a  clear  and  detailed  understanding  of  atomization  was  emphasized  in 
the  study  for  two  reasons:  (1)  atomization  provides  the  initial  conditions  for  subsequent 
combustion  processes  by  its  determinant  effect  on  drop  size  and  velocity;  and  (2)  the  periodic 
nature  of  primary  atomization  (ligament  shedding)  has  pronounced  similarities  to  combustion 
oscillations  in  rocket  engines  in  terms  of  both  frequency  range  and  the  frequency  dependency  on 
injector  operational  and  geometric  parameters. 

The  three  classical  cases  of  fully-developed  laminar  jet  flow,  fully-developed  turbulent  jet 
flow,  and  "plug"  jet  flow  were  studied.  These  cases  were  chosen  because  of  their  well- 
characterized  velocity  and  turbulence  intensity  profiles.  Absolute  plug  flow  conditions  were  not 
obtained,  but  were  approached  by  using  orifices  with  short  length-to-diameter  ratios  (Ud0~ 5). 
Length-to-diameter  ratios  in  impinging  jet  orifices  used  in  rocket  engines  are  typically  about 
three  to  five.  In  addition  to  changing  flow  condition  and  L/do,  other  test  parameters  (see  Table  1) 
included  impingement  angle,  20,  orifice  diameter,  and  free  jet  length  prior  to  impingement 

It  was  obvious  that  the  flow  condition  of  the  jet  before  impingement,  i.e.,  whether  it  was 
laminar  or  turbulent,  had  the  major  effect  on  atomization.  To  illustrate  the  importance  of  the 
initial  conditions  of  the  liquid  jet,  consider  the  instantaneous  images  of  the  sheets  formed  by 
laminar  and  turbulent  impinging  jets  shown  in  Fig.  2.  These  images  were  taken  under  quiescent 
conditions  and  at  atmospheric  pressure.  Although  the  jet  Reynolds  numbers  for  both  the  laminar 
and  turbulent  cases  are  similar,  the  resultant  sheets  have  very  different  and  distinct  character¬ 
istics.  In  the  laminar  case.  Fig.  2a,  small  ripples  on  the  surface  of  the  sheet  are  seen  near  the 
impingement  point,  and  after  some  distance  the  sheet  suddenly  disintegrates  into  drops.  In 
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(a)  (b) 

Fig.  2.  Instantaneous  images  of  sprays  formed  by  two  impinging  water  jets,  (a)  Resultant 
sheet  formed  by  two  laminar  impinging  jets  emanating  from  0.51  mm  inner  diameter, 

L/do=375  precision  bore  glass  tubes.  The  jet  velocity  was  13.1  m/s  and  the  full  impingement 
angle,  20,  was  60°.  (b)  Resultant  sheet  formed  by  two  turbulent  impinging  jets  emanating 
from  0.64  mm  inner  diameter,  L/do=80  precision  bore  glass  tubes.  The  jet  velocity  was 
12.2  m/s  and  the  full  impingement  angle,  20,  was  60°. 

many  of  the  images  of  laminar  impinging  jets,  incipient  breakup  occurred  at  mid-span  of  the 
sheet.  Drops  are  also  seen  shedding  steadily  off  the  edge  of  the  sheet.  In  the  turbulent  case.  Fig. 
2b,  the  waves  that  appear  at  impact  are  larger  than  in  the  laminar  case.  Downstream,  the  sheet 
disintegrates  into  ligaments,  with  incipient  breakup  consistently  occurring  at  the  edges  where  the 
sheet  is  thinnest  Periodicity  is  indicated  by  waves  on  the  liquid  sheet  and  by  the  spacing 
between  the  detached  ligaments.  Examination  of  the  detached  liquid  structures  reveals  irregular 
liquid  shapes  that  contract  into  roughly  cylindrical  ligaments.  The  cylindrical  ligaments  appear 
to  contract  further  into  irregularly-shaped  drops  that  will  eventually  take  a  nearly  spherical  shape. 
The  ligament-to-droplet  formation  process  at  atmospheric  pressure  is  most  likely  controlled  by 
surface  tension-driven  instabilities,  and  by  a  combined  surface  tension  -  aerodynamic  effect  at 
higher  ambient  pressures. 

The  turbulent  sheet  is  not  as  symmetric  as  the  laminar  sheet,  and,  by  comparing  upstream 
and  downstream  points  near  the  edge  of  the  turbulent  sheet  there  appears  to  be  relatively  large- 
scale  displacement  of  the  sheet  in  the  image  plane,  indicating  jet  unsteadiness  on  the  large  scale. 
Examination  of  opposite  edges  at  the  same  downstream  location  leads  to  the  same  conclusion 
regarding  asymmetry  and  large-scale  jet  unsteadiness.  This  phenomena  is  in  all  probability 
three-dimensional. 

From  the  results  presented  in  Fig.  2,  it  is  clear  that  the  phenomena  controlling 
atomization  for  impinging  jet  injectors  is  complicated  and  requires  both  theoretical  and 
experimental  study.  Thus,  emphasis  has  been  given  to  comparing  the  non-vaporizing 
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atmospheric  impinging  jet  studies  with  theoretical  models  as  well  as  extending  the  experimental 
work  to  higher  ambient  pressures  and  forced  oscillatory  conditions. 


Injector  Type 


Test  Parameter 

Glass  Tubes 

Twist-Prilled 

EDM 

Orifice  Diameter,  mm 

0.51, 0.64, 1.02, 1.45 

1.02 

0.51 

Impingment  Angle,  degrees 

40, 60, 80, 100, 180 

60 

60 

Length-to-Diameter  Ratio 

5, 35, 50,  80, 375 

10 

5 

Pre-impingement  Length,  mm 

2.5, 25, 35 

14 

2.5 

Laminar  Flow  Conditions: 

2800</te/<10000 

200cWe;<2300 

N/A 

N/A 

Turbulent  Flow  Conditions: 

4000</te/<30000 

3000<Refc20000 

4000</?ey<30000 

300<Wej<7000 

150<%<4000 

300<Wey<7000 

Ambient  Pressure 

1 

1,6,9 

1,5, 10 

Acoustic  Frequency,  Hz 

N/A 

670,1360,2020 

N/A 

Acoustic  Pressure,  kPa  (psi) 
(peak  -  peak  amplitude) 

N/A 

1.4 -6.2  (0.2 -0.9) 

N/A 
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2.1  Modeling  Results 

The  significant  task  of  developing  a  mechanistic  model  of  primary  atomization  was 
undertaken  as  part  of  this  contract  A  mechanistic  model  that  can  accurately  predict  the  effects 
of  injector  design  and  operation  on  breakup  length,  atomization  frequency,  and  drop  size  and 
velocity  distributions  is  critically  needed  by  combustor  designers.  Two  types  of  models  have 
been  studied  for  their  usefulness  for  predicting  impinging  jet  atomization:  (1)  a  linear  stability  - 
based  model  that  was  used  to  provide  predictions  of  breakup  length,  fastest  growing  wavelength, 
and  mean  drop  size  predictions;  and  (2)  a  finite  difference  model  that  was  used  to  study  the 
effects  of  jet  disturbances  and  jet  -  jet  interactions  on  incipient  wave  formation  on  the  liquid 
sheet 

2.1.1  Linear  Stability  Model 

The  atomization  model  that  has  been  used  to  date  by  most  workers  in  this  area4-5’6  is 
based  on  a  linear  stability  analysis  of  aerodynamically-induced  wave  growth  on  the  surface  of  a 
thinning,  viscous  liquid  sheet  This  model  was  used  to  predict  breakup  length,  the  periodic 
structure  of  breakup,  and  drop  size.  Details  of  our  implementation  of  the  model  as  well  as  a 
more  detailed  comparison  with  the  present  experimental  results  can  be  found  in  References  2  and 
3  and  Appendix  1. 

In  brief,  the  linear  stability  model  is  based  on  the  growth  of  infinitesimal  disturbances  due 
to  aerodynamic  stresses  on  the  liquid  sheet  surface  to  describe  the  disintegration  of  liquid  sheets. 
The  disturbance  on  the  sheet  surface,  rj,  is  given  by 

(1) 

1o 

where  7]0  is  the  initial  displacement  amplitude,  ft  is  the  growth  rate  and  t  is  time. 

Typically  the  growth  rate,  ft  ,is  calculated  for  a  spectrum  of  wavenumbers,  k.  The 
disturbance  wavenumber  corresponding  to  the  maximum  growth  rate,  ft  m,  controls  the  breakup 
process.  Both  sinuous  (antisymmetric)  and  dilatation  ol  (symmetric)  waves  can  grow;  however, 
previous  research  indicates  that  sinuous  waves  grow  faster  than  dilatational  waves,6  hence  only 
the  behavior  of  sinuous  disturbances  were  considered.  The  theory  does  not  predict  a  critical 
disturbance  amplitude  for  sheet  disintegration,  and  consequently,  an  empirical  relation  of  the 
following  form  is  typically  used:2 

h  H  ft  m 

fPi,mdt  =  /  7 f*  =  12  (2) 

0  0  us 
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where  tb  is  the  sheet  breakup  time,  xb  is  the  breakup  length,  and  Us  is  the  sheet  velocity.  The 
fastest  growing  wave  and  its  growth  rate  are  found  from  the  general  dispersion  equation: 


p2i,nd  + 


(of 

Res 


Pi,nd  + 


=  0 


(3) 


where; 


Wes 


=e£dh. 

a 


where  h  is  the  sheet  thickness  and  s  is  the  ratio  of  the  gas  density  to  the  liquid  density.  The  first 
term  in  Equation  (3)  is  associated  with  momentum  in  the  displacement  direction,  the  second  term 
rises  from  viscous  forces  while  the  two  parts  of  the  third  term  account  for  surface  tension  forces 
and  aerodynamic  forces,  respectively.  Drop  size  can  also  be  calculated  by  assuming  the  sinuous 
waves  that  are  formed  on  the  sheet  break  at  their  crests  and  troughs  and  subsequently  contract 
into  cylindrical  ligaments  with  diameter  i i.  Drops  are  formed  when  the  cylindrical  ligaments 
undergo  the  well-known  process  of  capillary  instability  and  pull  themselves  into  spheres  (drops) 
with  a  diameter,  do,  which,  under  the  conditions  of  the  present  study,  can  be  expressed  by 

dD~1.89dL  (4) 


Fig.  3.  Non-dimensional  lengths  as  a 
function  of  geometrically  scaled  jet  Weber 
number  as  predicted  by  aerodynamic  in¬ 
stability  model.  Lengths  normalized  by 
orifice  diameter. 


In  the  course  of  analyzing  results  from  the 
aerodynamic  stability  model,  a  non-dimensional 
scaling  parameter  based  on  the  jet  Weber  number, 
Wej-pjvjPdo/o,  and  the  half-impingement  angle,  0, 
was  identified:  Wej(l -cos  dp/sin3  6.  Use  of  this 
parameter  collapses  the  theoretical  dependence  of 
breakup  length,  fastest-growing  wavelength,  and 
drop  size  on  orifice  diameter,  impingement  angle, 
and  velocity  into  a  single  curve  dependent  on  the 
scaling  parameter.  Results  from  the  model  are 
shown  in  Fig.  3,  where  non-dimensional  lengths, 
normalized  by  the  orifice  diameter,  do,  are  plotted 
against  the  geometrically  scaled  jet  Weber  number. 
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The  formation  of  large  impact  waves  on  the  sheet  formed  by  impinging  turbulent  jets  is 
an  obvious  process  which  must  be  clearly  understood  before  an  accurate  mechanistic  atomization 
model  is  developed.  The  impact  waves  are  particularly  important  for  combustion  stability 
because  they  determine  the  temporal  nature  with  which  ligaments  are  shed  from  the  sheet,  and 
thus  control  the  temporal  nature  of  drop  formation  as  well.  The  wavelength  also  controls  the 
ligament  diameter  in  conjunction  with  the  thickness  of  the  sheet  at  its  breakup. 

The  observation  that  the  spacing  between  these  large  disturbances  is  primarily  dependent 
on  jet  diameter  and  apparently  independent  of  jet  velocity  suggests  that  the  source  of  the  impact 
waves  is  a  long- wave  instability  on  the  jet,  which  also  possesses  this  dependency.7  This  notion 
is  reinforced  by  the  observation  of  bulges  on  the  free  liquid  jet  prior  to  impingement  that  are 
similar  in  length  to  the  impact  waves.  A  simple  analytical  model  such  as  the  linear  stability 
model  described  in  Section  2.1.1  cannot  accurately  give  a  description  of  these  physical  processes. 
To  specifically  address  the  impact  wave  formation  process,  modeling  efforts  have  centered  on 
the  use  of  a  finite  difference  Navier-Stokes  code8  (RIPPLE)  to  investigate  the  effects  of  spatial 
and  temporal  jet  flow  oscillations  on  the  primary  breakup  of  the  impinging  jet  fan  into  ligaments. 
Directly  opposed  axisymmetric  jets  are  chosen  for  study  to  simplify  the  analysis. 

Representative  computational  results  showing  the  free  surface  contour  maps  of  two 
directly  opposed  impinging  jets  are  shown  in  Fig.  4.  Harmonic  disturbances  were  imposed  on 
the  jets  as  boundary  conditions.  The  cases  presented  include  imposed  sinusoidal  disturbances  of 
250, 500,  and  1000  s'1  on  each  jet  that  are  in  phase  with  each  other.  The  imposed  disturbance 
has  a  radial  profile  that  is  proportional  to  sin(nr/do),  with  the  maximum  axial  velocity  oscillation 
u'  at  the  jet  periphery  being  5%  of  the  mean  jet  velocity  Uj. 

The  computed  free  surface  contours  have  some  interesting  similarities  to  directly  opposed 
jets  observed  in  the  laboratory.  The  use  of  in-phase  jet  oscillations  result  in  the  antisymmetric 
sheet  disturbances  that  are  also  seen  in  the  laboratory.  The  computed  movement  of  disturbances 
on  the  liquid  sheets  has  been  analyzed  for  the  cases  with  500  and  1000  s*1  disturbances  imposed 
on  the  jet;  results  are  provided  in  Table  2  for  the  non-dimensionalized  disturbance  phase  speed, 
cAJ,  the  non-dimensionalized  disturbance  wavelength,  X/d0,  and  the  non-dimensionalized 
disturbance  frequency,  fshee/fjet.  It  is  seen  that  for  all  cases,  the  sheet  disturbance  frequency  is 
equal  to  the  frequency  of  the  imposed  oscillation  on  the  jet  Also,  the  non-dimensional 
disturbance  wavelength  for  the  500  s'1  imposed  oscillation  case  approaches  the  measured  non- 
dimensional  wavelength,  which  is  approximately  2.  This  model  has  some  promise  of  clarifying 
the  process  of  impact  wave  formation,  and  ongoing  work  is  focused  on  the  most  appropriate  way 
to  simulate  pre-impingement  jet  oscillations. 


Steady  inlet 
conditions 


250  Hz  Oscillation 


500  Hz  Oscillation  1000  Hz  Oscillation 


Fig.  4.  Computed  free  surface  contours  at  steady  inlet  conditions  and  with  harmonic  disturbances 
imposed  at  the  jet  inlet  Top  jet  enters  at  the  upper  left  hand  corner  with  a  mean  velocity  of  -10  m/s  and 
bottom  jet  enters  at  the  lower  left  hand  corner  with  a  mean  velocity  of  9.6  m/s.  The  jet  diameter  is  10 
mm  and  the  imposed  disturbances  are  5%  (maximum)  of  mean  velocity  with  sin(nr/2R)  profile. 


Table  2.  Summary  from  Finite-Difference  Computations  of  Unsteady  Impinging  Jets 


Frequency  of 

Lida 

f  sheet  jet 

Imposed  Oscillation.  Hz 

c/ll 

500 

0.60 

1.34 

1 

1000 

0.77 

0.78 

1 
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2.2  Experimental  Results 

A  series  of  detailed  measurements  and  photographic  and  electronic  images  were  made  to 
characterize  the  atomization  process  of  two  impinging  water  jets.  Both  laminar  and  turbulent  jets 
were  studied.  The  limiting  cases  of  fully-developed  laminar  and  fully-developed  turbulent  flow 
were  examined.  The  third  limiting  case  of  plug  flow  was  approached  by  using  orifices  with  an 
Ud0  of  5.  Some  attempts  were  made  to  study  impinging  jets  from  a  sharp-edged  orifice,  but 
were  cut  short  because  spiraling  instabilities  on  the  jets  could  not  be  avoided  in  spite  of  various 
means  taken  to  stabilize  the  flow.  Experimental  details  can  be  found  in  References  2  and  3. 
Table  1  summarizes  the  test  parameters. 

For  the  laminar  jets,  breakup  length  and  sheet  width  at  breakup  were  measured  from 
photographs  and  electronic  images  for  jets  emanating  from  0.51  mm  diameter  precision-bore 
glass  tubes  for  impingement  angles  of  40, 60,  80, 100,  and  180°.  Tests  were  done  at  open  air 
conditions.  The  orifice  length  was  190  mm  ( Udo  =  375),  resulting  in  fully-developed  laminar 
flow  at  the  orifice  exit  The  free  jet  length  prior  to  impingement  was  10  mm.  Flow  Reynolds 
numbers  from  2800  to  9100  corresponding  to  jet  velocities  from  5.5  to  18  m/s  were  achieved; 
laminar  flow  at  the  higher  velocities  was  possible  due  to  a  contoured  inlet  and  the  smooth  glass 
tubes.  Weber  numbers,  Wej,  varied  from  200  to  2200. 

Because  the  practical  interest  is  in  turbulent  impinging  jets,  their  experimental 
characterization  was  much  more  extensive  than  it  was  for  laminar  impinging  jets.  The 
experimental  parameters  were  extended  to  include  effects  of  orifice  diameter,  orifice  length,  free 
jet  impingement  length,  impingement  angle,  and  jet  velocity.  Furthermore,  to  evaluate  the 
effects  of  practical  fabrication  techniques,  impinging  jets  from  twist-drilled  and  electro-discharge 
machined  (EDM)  orifices  were  studied  in  addition  to  the  precision-bore  glass  tubes.  The  twist- 
drilled  orifices  were  1.0  mm  in  diameter,  10  mm  long,  and  produced  free  jets  that  were  14  mm 
long  prior  to  impingement  The  EDM  injector  was  0.51  mm  in  diameter  and  was  2.5  mm  long 
with  a  free  jet  length  of  2.5  mm  before  impingement  For  both,  the  impingement  angle  was  60°. 
Measurements  included  breakup  length,  sheet  width  at  breakup,  the  distance  between  apparent 
disturbances  on  the  liquid  sheet  surface,  and  drop  size  and  drop  velocity  distributions.  For  the 
twist-drilled  and  EDM  injectors,  the  above  measurements  were  also  made  in  closed  chambers  at 
elevated  ambient  pressures,  and  under  forced  oscillatory  ambient  conditions  for  the  twist-drilled 
injector. 
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221 _ Atmospheric  Pressure  Atomization  Studies  of  Impinging  Jets 

from  Precision-Bore  Glass  Tubes 

Extensive  experimental  characterization  was  made  at  atmospheric  conditions  with  the 
precision-bore  glass  tube  injectors  for  both  laminar  and  turbulent  jet  conditions  to  provide  a 
sound  basis  for  later  comparisons  with  data  from  more  practical  injector  types  and  for  data  taken 
under  high-pressure,  acoustic,  and  combusting  flow  conditions.  Most  of  the  glass  tube  results  are 
detailed  in  References  2  and  3  and  in  Appendix  1.  A  brief  summary  of  the  breakup  length,  drop 
size  distribution,  and  disturbance  wavelength  measurements  are  provided  here. 

Length  measurements  of  the  intact  sheet  and  of  periodic  structures  were  made  from 
images  such  as  those  shown  in  Fig.  2.  Non-dimensional  breakup  length,  xi/do,  is  plotted  against 
jet  Weber  number  in  Fig.  5  for  laminar  and  turbulent  impinging  jet  cases.  Clear  differences 
between  the  two  cases  can  be  seen.  The  non-dimensional  breakup  length  of  the  turbulent 
impinging  jets  appears  to  have  some  dependence  on  impingement  angle,  but  appears  to  be 
relatively  independent  of  jet  Weber  number.  There  is  a  strong  dependence  on  jet  Weber  number 
for  the  laminar  impinging  jet  case. 


Fig.  5a.  Non-dimensional  breakup 
length  of  laminar  impinging  jets  from 
precision  bore  glass  tubes  as  a  function 
of  jet  Weber  number. 


Fig.  5b.  Non-dimensional  breakup 
length  of  turbulent  impinging  jets  from 
precision  bore  glass  tubes  as  a  function 
of  jet  Weber  number. 


Comparisons  of  the  measured  non-dimensional  breakup  length  with  predictions  from  the 
aerodynamic  instability  model  described  in  Section  2.1.1  are  shown  in  Fig.  6.  Again  the 
differences  between  the  turbulent  impinging  jet  case  and  the  laminar  impinging  jet  case  are 
evident.  Laminar  impinging  jet  breakup  may  be  modeled  relatively  well  by  the  aerodynamic 
instability  model,  but  again,  the  turbulent  impinging  jet  breakup  shows  little  dependence  on  the 


geometrically  scaled  jet  Weber  number.  It  appears  that  turbulent  impinging  jet  breakup,  which  is 
clearly  dependent  on  orifice  diameter  and  impingement  angle,  is  related  to  the  sheet  thickness. 


which  also  is  determined  by  orifice  diameter  and  impingement  angle . 
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The  measured  distance  normalized  by 
orifice  diameter  of  the  separation  between  periodic 
structures  for  turbulent  impinging  jets  is  shown  in 
Fig.  7  as  a  function  of  jet  velocity.  The 
"wavelength"  appears  to  be  primarily  dependent  on 
orifice  diameter  and  independent  of  jet  velocity. 
The  observation  of  independence  from  jet  velocity 
is  contrary  to  the  aerodynamic  instability  model 
predictions.  The  distance  between  detached 
ligaments  is  approximately  constant  at  around  4 
orifice  diameters,  and  the  distance  between  surface 
wave  structures  is  constant  at  around  2  orifice 
diameters.  There  was  a  rather  large  spread  in  the 
data,  with  a  standard  deviation  of  +/-  35%. 

Of  course,  the  combustion  process  must 
also  be  considered  in  an  analysis  of  combustion 
instability.  In  liquid  rocket  engine  combustion, 
vaporization  is  the  rate-limiting  step.  The 
vaporization  process  is  controlled  in  large  by  the 
drop  size.  Accurate  measurements  of  the  drop  size 
distribution  are  necessary  to  develop  an  accurate 
understanding  of  the  problem  of  combustion 
instability.  An  argon-ion  based,  two-component 
Phase  Doppler  Particle  Analyzer  (PDPA)  was 
used  for  making  drop  size  and  velocity 
measurements.  A  description  of  the  theoretical 
and  operating  principles  of  the  PDPA  are  given 
elsewhere.3  Non-dimensional  arithmetic  mean 
drop  diameters  ( Pu/do )  are  presented  in  Fig.  8 
for  the  turbulent  impinging  jet  case.  A 
comparison  of  measured  non-dimensional 
arithmetic  mean  drop  diameter  with  the 
monodisperse  drop  size  predicted  by  the 
aerodynamic  instability  model  is  also  shown. 

The  trend  of  the  model  predictions  match  the 


Fig.  6.  Non-dimensional  breakup  length 
as  a  function  of  geometrically  scaled  jet 
Weber  number.  Comparison  of  measurements 
with  aerodynamic  instability  model  predictions. 


Fig.  7.  Non-dimensional  distance 
measured  between  adjacent  surface  wave 
structures  and  between  adjacent  ligaments. 


•in^ 


Fig.  8.  Non-dimensional  drop  diameter  as  a 
function  of  geometrically  scaled  jet  Weber 
number.  Comparison  of  PDPA  Measurements 
with  predictions  from  aerodynamic  instability 
model. 
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experimental  trend  quite  well,  and  the  measurements  appear  to  collapse  into  a  single  curve  when 
plotted  against  the  geometrically  scaled  jet  Weber  number  as  the  model  predicts.  Adjustment  to 
the  empirical  breakup  constant  could  yield  a  very  good  match.  However,  based  on  the 
significant  disagreements  between  prediction  and  measurement  for  wavelength  and  breakup 
length,  the  mean  drop  size  agreement  is  thought  to  be  coincidental. 

An  empirical  correlation  for  normalized  drop  size  as  a  function  of  the  geometrically 
scaled  jet  Weber  number,  Wejfld),  was  obtained  with  a  correlation  coefficient  of  0.964: 

~~  =2.217  ■  {Wejf(d)}~°  354  (5) 

“o 

2.22  Effects  of  Ambient  Conditions  and  Ipjector  Type 

Experiments  using  precision-bore  glass  tubes  provided  fully-developed  laminar  and 
turbulent  jets  with  well-characterized  velocity  and  turbulence  intensity  profiles.  In  practical 
combustors,  however,  injection  orifices  and  impingement  lengths  are  short,  typically  less  than  5, 
and  ambient  pressures  are  high  and  often  unsteady.  To  evaluate  these  effects,  measurements 
were  also  taken  using  twist-drilled  and  EDM  injectors  and  at  high  pressure  and  oscillatory 
ambient  conditions  in  a  transparent  acoustic  chamber.  As  with  the  quiescent,  atmospheric 
pressure  case,  measurements  made  within  the  confines  of  the  acoustic  chamber  included  sheet 
breakup  length  and  drop  size  distribution. 

Instantaneous  images  of  the  EDM  injector's  spray  at  1, 5,  and  10  atm  ambient  pressure 
are  shown  in  Fig.  9.  These  images  clearly  show  that  faster  ligament  breakup  and  drop  dispersion 
result  from  increases  of  ambient  pressure.  A  plot  of  the  non-dimensional  breakup  length,  Xf/d0, 
as  a  function  of  Wej  (=pjVJ^d(/a)  made  from  such  images  is  shown  in  Fig.  10.  Each  symbol  in 
Fig.  10  represents  an  average  of  17  individual  breakup  length  measurements,  while  the  bars 
represent  the  plus/minus  standard  deviation  of  those  measurements.  Reductions  in  breakup 
length  are  seen  with  increasing  ambient  pressure.  It  is  also  seen  that,  as  for  the  case  of  glass 
tubes  at  atmospheric  ambient  conditions,  increases  in  Wej  apparently  have  little  effect  on  the 
non-dimensional  breakup  length. 

One  of  the  major  differences  between  the  glass  tubes  and  the  EDM  injector  is  die  orifice 
length-to-diameter  ratio,  L/d0.  The  glass  tube  L/d0  is  80  while  that  of  the  EDM  injector  is  about 
5.  In  spite  of  the  large  difference  in  Ud0,  there  is  no  apparent  difference  between  the  measured 
non-dimensional  breakup  lengths  for  the  glass  tube  injector  and  the  EDM  injector.  This 
observation  that  the  length-to-diameter  ratio  does  not  appreciably  affect  the  impinging  jet  spray 
breakup  length  was  also  noted  by  Anderson,  et  al.2  Another  difference  between  the  two  injectors 
was  the  length  of  the  free  jet  prior  to  impingement:  for  the  glass  tube  it  was  25.4  mm,  and  for  the 
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EDM  injector  it  was  2.5  mm.  Anderson,  et  al.2  also  noted  that  free  jet  length  had  no  apparent 
effect  on  spray  formation  processes  of  impinging  turbulent  liquid  jets. 


(a)  (b)  (c) 

Fig.  9.  Instantaneous  images  of  the  spray  formed  by  an  electro-discharge  machined  (EDM)  impinging 
jet  injector.  Orifice  diameter  was  0.51  mm,  and  impingment  angle  was  60°.  Ambient  pressure  was 
(a)  1  atm;  (b)  5.1  atm;  (c)  10.5  atm. 


Fig.  10.  Non-dimensional  breakup  length 
as  a  function  of  the  jet  Weber  number, 
ambient  pressure,  and  injector  type. 


Fig.  11.  Non-dimensional  surface 
disturbance  wavelength  as  a  function 
of  the  jet  Weber  number,  ambient  pressure, 
and  injector  type. 
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The  effects  of  ambient  pressure  and  injector  type  on  non-dimensional  liquid  sheet 
disturbance  wavelength  are  shown  in  Fig.  11.  A  small  reduction  in  disturbance  wavelength 
appears  with  the  EDM  injector,  which  has  a  substantially  smaller  orifice  length  and  a 
substantially  smaller  free  jet  length  than  does  the  glass  tube.  However,  as  with  the  glass  tubes, 
the  non-dimensional  disturbance  wavelength  appears  to  be  essentially  independent  of  velocity,  or 
of  jet  Weber  number  as  plotted  in  Fig.  1 1.  The  non-dimensional  disturbance  wavelength  also 
appears  to  be  independent  of  ambient  pressure,  at  least  to  10  atm,  which  further  confirms  the 
inadequacy  of  the  aerodynamic  instability  model  at  these  conditions. 

Drop  size  data  for  the  EDM  injector  at  atmospheric  and  higher  ambient  pressures  are 
compared  to  glass  tube  injector  drop  size  data  at  atmospheric  ambient  pressure  in  Fig.  12.  A 
small  reduction  in  non-dimensional  drop  size  is  observed  for  the  EDM  injector  vs  the  glass  tube 
injector  at  atmospheric  ambient  pressure,  and  another  small  reduction  is  observed  when  the  EDM 
injector  is  operated  at  4.6  atm  ambient  pressure. 


Fig.  12.  Non-dimensional  arithmetic  drop  size 
measured  at  the  impinging  jet  spray  centerline 
as  a  function  of  the  geometrically-scaled  jet 
Weber  number  and  injector  type.  Drop  size 
measurements  for  the  EDM  injector  were  taken 
at  1  and  4.6  atmospheres. 


Together,  the  data  presented  in  Figs.  10, 1 1,  and  12  can  be  used  to  conjecture  a  partial 
mechanism  for  the  atomization  of  liquid  sheets  formed  by  impinging  jets.  The  impact  waves 
travel  out  from  the  impingement  point  down  the  thinning  sheet  with  an  undetermined  phase 
speed.  Incipient  ligament  formation  occurs  due  to  combined  aerodynamic  and  tensile  forces, 
which  accounts  for  the  effect  of  ambient  pressure  on  breakup  length;  the  impact  wave  formation 
process  in  not  considered  to  be  affected  by  aerodynamic  pressure.  The  ligament  then  fragments, 
also  under  combined  aerodynamic  and  tensile  forces.  It  is  not  clear  why  ligament  spacing  is 
twice  the  spacing  between  adjacent  surface  disturbances.  Although  the  effect  of  increased 
ambient  pressure  on  the  ligament  fragmentation  process  is  to  produce  smaller  drops,  it  also 
results  in  a  shortened  breakup  length,  and  thus  larger  ligaments  from  the  thicker  sheet  because 
the  impact  disturbance  wavelength  does  not  change.  The  net  effect  is  a  small  change  in  drop  size 
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due  to  ambient  pressure  effects.  This  phenomenological  model  is  completely  consistent  with  the 
data  shown  in  Figs.  10, 11,  and  12. 

The  source  of  the  impact  waves  remains  a  key  question.  Because  their  characteristic 
length  is  apparently  independent  of  jet  velocity,  impingement  angle,  or  ambient  pressure,  and 
linearly  dependent  on  jet  diameter,  it  may  be  postulated  that  the  large  impact  waves  seen  on  the 
sheet  within  a  few  diameters  downstream  of  the  impingement  point  have  their  origin  in  the  jet 
upstream  of  the  impingement  point  Instabilities  seen  on  the  free  jet  appear  to  have  a 
characteristic  length  that  is  close  to  twice  the  jet  diameter,  the  same  as  the  surface  impact  waves. 
It  should  be  noted  that  helical  instabilities  on  high  velocity  liquid  jets  are  predicted  to  have 
precisely  these  characteristics,  i.e.,  the  fastest  growing  disturbance  has  a  wavelength 
approximately  equal  to  twice  the  jet  diameter,  and  is  independent  of  jet  velocity.7  It  is  not  clear 
whether  there  is  an  interaction  mechanism  that  occurs  between  the  jets.  An  alternative 
explanation  for  the  source  of  the  impact  waves  is  that  they  are  generated  from  oscillations 
inherent  in  stagnation  point  flow. 

To  examine  the  effects  of  an  oscillating  pressure  field,  an  acoustic  chamber  (Fig.  13)  was 
designed  to  excite  the  first  (1W),  second  (2W)  and  third  (3W)  resonant  modes  at  frequencies  on 
the  order  of  several  thousand  hertz  under  ambient  temperature  conditions  for  air  and  helium 
environments.  The  interior  dimensions  of  the  chamber  are  254  mm  in  width,  305  mm  in  height 
and  102  mm  in  depth.  Two  large  Plexiglas  faces  allow  visual  access  into  the  chamber,  while 
steel  grids  brace  the  Plexiglas  for  pressurized  studies.  The  maximum  chamber  pressure  is 
approximately  1.03  MPa  (150  psia).  These  modes  and  frequencies  were  chosen  because  of  their 
similarity  to  typical  rocket  instabilities.  This  frequency  range  also  allows  for  tests  in  stable  and 
unstable  regions  according  to  the  Hewitt  stability  correlation  for  realistic  injector  stability 
parameter  values. 


Fig.  1 3.  An  image  of  the  transparent  acoustic 
chamber.  The  interior  dimensions  of  the 
chamber  are  254  mm  in  width,  305  mm  in 
height,  and  102  mm  in  depth.  An  acoustic 
compression  driver  is  encased  in  a  steel 
cylindrical  pressure  vessel  attached  to  the 
right  side-wall  of  the  chamber.  The 
impinging  jet  injector  is  inserted  into  the 
top  plate  of  the  chamber. 


An  acoustic  driver  with  its  highest  efficiency  in  the  thousand  hertz  range  was  attached  to 
the  chamber  right  side-wall  and  encased  in  a  steel  pressure  vessel  to  keep  equal  pressure  across 
the  acoustic  driver  diaphragm  during  pressurized  experiments.  The  speaker  was  driven  at  a 
particular  frequency  and  a  microphone  was  used  to  measure  the  oscillatory  pressure  field  within 
the  chamber.  The  twist-drilled  or  EDM  impinging  element  fit  into  the  top  plate  of  the  chamber. 
Water  flowed  through  a  41.3  mm  inner  diameter  stilling  chamber  fitted  with  a  sintered  brass  disc 
before  entering  the  injector  orifices. 

To  characterize  the  dimensionality  of  the  applied  acoustic  field,  extensive  acoustic 
measurements  were  taken  at  the  three  resonant  modes.  The  excited  1W,  2W  and  3W  modes, 
which  had  frequencies  of  668  Hz,  1362  Hz  and  2016  Hz,  respectively,  were  primarily  one¬ 
dimensional  with  minor  variations  in  wave  amplitude  from  the  top  to  bottom  of  the  chamber. 

Fig.  14  shows  the  measured  pressure  amplitude  as  a  function  of  position  across  the  chamber 
width  at  atmospheric  pressure.  The  shape  of  the  waveforms  are  very  similar  to  the  expected  first 
through  third  standing  wave  modes.  All  three  waveforms  were  generated  with  the  same  input 
voltage  and  amplifier  gain  to  the  speaker.  The  peak  pressure  amplitude  for  the  2W  mode  is 
approximately  3.0  kPa  (0.44  psi)  at  the  center  of  the  chamber,  and  is  greater  than  the  other  two 
modes.  The  maximum  observed  pressure  amplitude  is  approximately  1.5  kPa  (0.22  psi)  for  the 
1 W  mode  and  0.70  kPa  (0. 10  psi)  for  the  3W  mode. 


Fig.  14.  The  measured  acoustic  pressure 
amplitude  plotted  as  a  function  of  position 
along  the  chamber  width  for  the  1W,  2W, 
and  3W  resonant  modes  at  a  chamber 
pressure  of  one  atmosphere. 


The  effects  on  spray  formation  of  oscillating  velocity  and  of  oscillating  pressure  can  both 
be  evaluated  at  the  chamber  centerline  where  the  injector  is  located  by  exciting  the  1 W  and  2W 
mode,  respectively.  Since  the  intended  resonant  acoustic  mode  structures  are  satisfactory,  the 
effects  of  the  applied  acoustic  field  on  the  impinging  jet  spray  and  their  relation  to  the  Hewitt 
stability  correlation  can  be  pursued. 

Curtains  of  air  across  portions  of  the  chamber's  Plexiglas  walls  were  used  to  keep  water 
wetting  to  a  minimum.  Breakup  length  measurements  could  not  be  obtained  at  higher  Weber 
numbers  and  at  higher  chamber  pressures  because  of  an  inability  to  acquire  high  contrast  images 
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due  to  the  nature  of  the  turbulent  spray  as  well  as  lighting  problems  caused  by  excessive  window 
wetting.  Improvements  in  window  cleaning  and  lighting  practices  are  needed  to  make 
measurements  at  higher  Weber  numbers  at  pressurized  conditions. 

Drop  size  and  velocity  measurements  were  made  using  the  PDPA  along  the  twist-drilled 
injector  spray  centerline  which  coincides  with  the  chamber  centerline  at  a  distance  of  25  mm 
downstream  of  the  impingement  point  Axial  and  transverse  (across  the  chamber  width)  drop 
velocities  were  measured.  The  injection  velocity  of  the  water  was  9.0  m/s.  Five  thousand  five 
hundred  valid  drop  size  and  velocity  samples  were  obtained  at  this  point  for  each  test  condition. 
The  effects  of  an  applied  acoustic  field  were  investigated  by  subjecting  the  impinging  jet  spray  to 
the  1W  mode  at  668  Hz.  A  pressure  node  (oscillating  velocity)  existed  at  the  chamber  centerline 
where  the  measurements  were  made.  The  applied  1W  acoustic  field  did  not  have  an  appreciable 
effect  on  drop  size  or  drop  velocity  at  the  prescribed  measurement  location  for  a  chamber 
pressure  of  one  atm,  but  did,  however,  noticeably  affect  the  drop  formation  process  at  a  higher 
chamber  pressure  of  six  atm,  i.e.,  higher  ambient  densities  enhanced  the  effects  of  the  applied 
acoustic  field.  A  smaller  arithmetic  mean  drop  size,  approximately  230  pm,  was  measured  for 
the  acoustically  perturbed  case,  as  compared  to  the  quiescent  case,  where  D]g=  250  pm.  A  plot 
of  the  drop  size  number  distribution,  fiD),  versus  measured  drop  diameter,  D,  for  acoustically 
and  non-acoustically  perturbed  cases  is  illustrated  in  Fig.  15.  Evident  from  Fig.  15  is  the  higher 
drop  number  distribution  between  diameters  of  100  to  200  pm  for  the  acoustically  perturbed 
case.  This  contributes  to  the  observed  lower  arithmetic  mean  diameter.  Results  for  the  3W  mode 
were  similar  to  the  1W  mode  results,  whereas  the  results  from  an  applied  2W  mode  (pressure 
anti-node  at  spray  centerline)  indicated  that  oscillating  pressure  had  no  apparent  effects. 


Fig.  15.  Drop  size  number  distribution, 
flJD),  measured  by  the  PDPA  at  25  mm 
downstream  of  the  impingement  point 
without  and  with  an  applied  acoustic  field 
(1 W  mode).  Ambient  pressure  is 
approximately  6  atmospheres. 
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2.2.3  Combustine-Flow  Chamber  Design 

Definitive  information  regarding  combustion  instability  mechanisms  can  only  be 
obtained  under  high-pressure,  combusting-flow  conditions.  Furthermore,  it  is  very  desirable  to 
be  able  to  study  combustion  phenomena  under  conditions  of  oscillating  pressure  fields  similar  to 
those  that  actually  occur  during  a  combustion  instability  event.  A  high-pressure,  optically- 
accessible  combustor  in  which  oscillating  pressure  fields  can  be  generated  has  been  designed  to 
enable  experimental  studies  under  simulated  combustion  instability  conditions  at  the  Cryogenic 
Combustion  Laboratory.  Liquid  hydrocarbon  fuel  test  capability  for  impinging  jet  combusting- 
flow  tests  was  added  to  the  Cryogenic  Combustion  Laboratory.  A  schematic  of  the  chamber  is 
provided  in  Fig.  16.  The  short  orifice  (L/do=5),  short  impingement  length  (Limp/d0=5)  EDM 
injector  unit  that  matches  practical  injector  designs  and  that  has  undergone  cold-flow 
characterization  will  be  used  in  the  combusting-flow  tests. 


•  Nominal  operating  conditions 

-  Pressure  =  250  psia 

-  Gas  temperature  =  1500  K 

-  Resonant  frequency  =  2000  H 


•  Chamber  dimensions  (m) 
-  0.2  w;  0.25  h;  0.1  d 


•  Measurements 

-  High  speed  cinematography 

-  Photography 

-  Phase-locked  emission 

-  Two  component  PDPA 

-  HF  pressure 


•  Injector  location 
-  Pressure  Node/ 

Antinode  •  Side-mounted  drivers 


NOZZLE 


Fig.  16.  Schematic  of  high-pressure,  combusting-flow,  optically-accessible  combustor  with  features. 
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2.3  Discussion  of  Potential  Instability  Mechanisms 

A  number  of  theories  and  observations  have  been  presented  that  have  attempted  to 
explain  the  mechanism  of  combustion  instability.  A  thorough  review  of  past  work  done  in 
combustion  instability  was  done  under  this  contract,  and  is  summarized  in  Reference  1.  Because 
the  general  objective  of  this  work  was  to  develop  a  rational  physical  basis  for  the  Hewitt  Stability 
Correlation,  the  experimental  results  obtained  to  date  will  now  be  examined  in  the  context  of  the 
stability  correlation.  Specifically,  the  dependence  of  breakup  length,  drop  size  distribution,  and 
disturbance  wavelength  (frequency)  on  d(/Uj  will  be  examined.  It  is  important  to  keep  in  mind 
that  the  discussion  to  follow  pertains  to  cold-flow  atomization  data,  and  that,  ideally,  the 
functional  dependence  of  the  impinging  jet  atomization  characteristics  should  be  evaluated  under 
actual  rocket  conditions.  To  aid  understanding  the  discussion  to  follow,  it  is  useful  to  remember 
than  an  increase  in  the  stability  parameter,  dJUj,  leads  to  an  increase  in  stability. 

The  functional  dependence  of  the  measured  breakup  length  for  the  glass  tube  injector  on 
the  impingement  angle  and  dg/Uj  is  shown  in  Fig.  17  for  an  orifice  diameter  of  0.64  mm.  For  all 
impingement  angles,  the  nondimensional  quantity,  x//d0,  decreases  with  increasing  dJUj . 
Similar  behavior  is  noted  for  different  orifice  diameters,  as  displayed  in  Fig.  18;  however,  with 
an  increase  in  orifice  diameter,  the  nondimensional  breakup  length  is  less  sensitive  to  variations 
in  d(/Uj .  The  trend  of  the  measured  breakup  length  as  a  function  of  chamber  pressure  and 
injector  stability  parameter  is  shown  in  Fig.  19  for  the  EDM  injector  (d0= 0.51  mm).  Under 


Fig.  17.  Non-dimensional  breakup  length 
of  liquid  sheet  formed  by  impinging  jets  from 
glass  tubes  as  a  function  of  stability 
parameter  and  impingement  angle. 


Fig.  18.  Non-dimensional  breakup  length 
of  liquid  sheet  formed  by  impinging  jets  from 
glass  tubes  as  a  function  of  stability 
parameter  and  orifice  diameter. 
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atmospheric  pressure  (,PC=101  kPa),  the  nondimensional  breakup  length  initially  increases  then 
decreases  with  increasing  d^Uj .  With  increasing  chamber  pressure,  the  nondimensional 
breakup  length  seems  less  sensitive  to  changes  in  d(/Uj .  By  itself,  the  observation  that  breakup 
length  decreases  with  increasing  df/Uj  is  not  consistent  with  the  notion  that  engine  stability 
increases  with  an  increase  in  the  injector  stability  parameter.  A  reduction  in  sheet  breakup  length 
most  likely  moves  the  combustion  zone  closer  to  the  injector  face  plate,  thus  increasing  the 
probability  of  the  initiation  of  instability  phenomena.  It  should  be  noted  however  that  the  sheet 
breakup  length  stays  approximately  constant  with  increasing  d^Uj  at  higher  ambient  pressures 
as  shown  in  Fig.  19. 


Fig.  19.  Non-dimensional  breakup  length 
of  liquid  sheet  formed  by  impinging  jets 
from  EDM  injector  as  a  function  of 
stability  parameter  and  ambient  pressure. 


The  mean  drop  size  and  the  distribution  of  drop  size  may  be  critical  in  terms  of  whether 
instability  will  ensue.  Figure  20  shows  the  non-dimensional  arithmetic  mean  diameter,  D1(/d0, 
along  the  spray  centerline  at  an  axial  location  of  41  mm  from  the  impingement  point,  plotted  as  a 
function  of  dJUj  and  29.  Note  that  the  orifice  diameter  is  constant,  so  changes  in  dJUj  result 
only  from  the  variation  in  jet  velocity.  An  increase  in  mean  drop  size  is  noted  with  increasing 
d(/Uj ,  thus  leading  to  the  conclusion  that  bigger  drops  have  a  stabilizing  effect,  which  has  been 
widely  suggested,  but  never  substantiated.1  Similar  mean  drop  size  trends  with  dJUj  are 
observed  for  different  orifice  diameters  and  at  different  locations  in  the  spray  field.  It  should  be 
noted  that  the  same  mean  drop  size  will  not  necessarily  be  obtained  for  any  given  value  of  df/Uj . 

Another  important  consideration  is  the  degree  of  polydispersity  of  the  drop  size 
distribution.  It  is  reasonable  to  expect  that  stability  will  be  enhanced  if  there  is  a  wide 
distribution  of  drop  sizes  because  any  present  effects  of  resonant  burning  can  be  essentially 
neutralized  by  different-sized  drops  that  release  most  of  their  chemical  energy  out  of  phase  with 
the  drops  that  are  burning  in  resonance  with  pressure  oscillations.  Conversely,  in  the  case  of  a 
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relatively  monodisperse  drop  distribution,  the  drops  will  bum  similarly,  i.e.,  have  a  similar  drop 
lifetime.  The  number  distribution,  ./(dp),  plotted  as  a  function  of  drop  size,  dp,  and  the  injector 
stability  parameter  is  shown  in  Fig.  21.  The  measured  distributions  were  taken  along  the  spray 
centerline  (x=0)  at  an  axial  location  (z)  of  41  mm  for  an  injector  orifice  diameter  of  0.64  mm  and 
an  impingement  angle  of  60°.  The  change  in  d(/Uj  was  obtained  here  by  changing  jet  velocity 
while  maintaining  orifice  diameter.  Notice  that  an  increase  in  the  injector  stability  parameter 
leads  to  broader  distributions  with  larger  mean  values.  Similar  distribution  shapes  are  observed 
for  different  radial  positions  and  impingement  angles  for  a  given  orifice  diameter. 

Changes  in  the  stability  parameter  can  also  be  realized  by  changing  d0.  Eight  number 
distributions  as  a  function  of  dp  and  d</Uj  are  shown  in  Fig.  22  for  two  different  injectors.  The 
hollow  symbols  represent  data  from  the  0.64  mm  diameter  glass  tube  injector  and  the  solid 
symbols  represent  data  from  the  1.45  mm  diameter  glass  tube  injector.  All  the  data  was  acquired 
along  the  spray  centerline  and  at  a  similar  non-dimensional  axial  location.  Once  again,  an 
increase  in  df/Uj  -  this  time  obtained  via  a  combination  of  changes  in  do  and  Uj  -  leads  to 
broader  and  flatter  distributions  with  a  corresponding  increase  in  mean  drop  size. 

Finally,  it  is  clear  that  periodic  drop  formation  could  lead  to  periodic  burning  and 
subsequent  instability  if  the  frequency  with  which  the  drops  entered  the  combustion  zone 
matched  in  some  fashion  some  resonant  frequency  of  the  chamber.  Disturbance  wavelength  data 
was  converted  to  atomization  frequency  simply  by  dividing  the  jet  velocity  by  the  measured 
wavelength;  the  assumption  that  the  disturbance  velocity  is  equal  to  the  sheet  velocity  needs  to 
be  substantiated.  Earlier  studies2  showed  that  measured  drop  velocities,  anu  therefore  the  sheet 
velocity,  were  nearly  equal  to  the  jet  velocity.  Atomization  frequency  data  is  presented  in  Fig. 

23  based  on  ligament  spacing.  The  trend  of  atomization  frequencies  are  quite  similar  to  the 
maximum  possible  combustion  instability  frequency  given  by  the  Hewitt  Stability  Correlation, 
also  shown  in  Fig.  23.  The  similar  dependencies  of  maximum  possible  instability  frequency  and 
atomization  frequency  on  d(/Uj  is  significant  in  terms  of  primary  atomization  being  a  key 
process  in  combustion  instability.  The  calculated  atomization  frequency  is  about  two  times 
greater  than  the  maximum  combustion  instability  frequency  as  predicted  by  the  stability 
correlation,  however.  The  effects  of  periodic  atomization  must  be  more  carefully  analyzed,  and 
the  atomization  frequency  must  still  be  determined  under  high  pressure,  combusting-flow 
conditions.  There  must  also  be  an  accounting  for  the  effects  of  drop  size  distribution  as  shown  in 
Figs.  21  and  22. 
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Fig.  20.  Non-dimensional  arithmetic  mean  drop 
size  of  spray  formed  by  turbulent  impinging  jets 
from  glass  tubes  as  a  function  of  stability 
parameter  and  impingement  angle. 


Fig.  21.  Number  distribution  of  drop  sizes  in 
spray  formed  by  turbulent  impinging  jets  from 
glass  tubes  as  a  function  of  die  stability 
parameter. 
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Fig.  22.  Number  distribution  of  drop  sizes  in 
spray  formed  by  turbulent  impinging  jets  from 
glass  tubes  as  a  function  of  the  stability 
parameter.  Data  from  two  different 
injectors  are  shown. 


Fig.  23.  Atomization  frequency  of  spray 
formed  by  turbulent  impinging  jets  from 
glass  tubes  as  a  function  of  stability 
parameter.  Frequency  is  assumed  to  be  equal 
to  the  ratio  of  sheet  speed  to  distance  between 
adjacent  detached  liquid  ligaments.  Also 
shown  is  maximum  instability  frequency 
predicted  by  the  stability  correlation. 
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3.0  SUMMARY 

The  most  complete  and  systematic  study  of  the  atomization  of  impinging  liquid  jets  to 
date  has  been  completed  under  this  contract  Effects  of  jet  flow  condition,  orifice  diameter, 
impingement  angle,  pre-impingment  length,  fabrication  procedure,  and  jet  velocity  at  steady  and 
oscillating,  and  atmospheric  -  and  high  -  pressure  conditions  have  been  investigated.  Results  of 
these  experimental  studies  have  been  compared  to  current  theories  in  terms  of  sheet  breakup 
length,  wave  structure,  and  drop  size.  Experiments  contrasting  laminar  and  turbulent  jet 
conditions  clearly  demonstrate  that  the  jet  conditions  have  a  dramatic  effect  on  the  atomization 
process.  Specifically,  the  measured  breakup  lengths  for  laminar  impinging  jets  are  longer  and 
displayed  different  trends  as  a  function  of  jet  velocity  and  impingement  angle  than  those  for 
turbulent  impinging  jets.  In  addition,  the  measured  spray  width  for  the  turbulent  impinging  jet 
case  remained  approximately  constant  as  a  function  of  jet  velocity  and  impingement  angle  for  a 
fixed  orifice  diameter,  unlike  the  laminar  impinging  jet  case.  The  present  results  are  generally  in 
agreement  with  earlier  studies  at  atmospheric  pressure  by  Dombrowski  and  Hooper,10  however, 
there  are  important  discrepancies  with  their  speculation  that  these  differences  are  due  to  velocity 
profile  differences  between  fully-developed  laminar  flow  and  fully-developed  turbulent  flow; 
implicit  in  their  reasoning  was  the  assumption  that  aerodynamic  forces  played  a  large  role  in  the 
formation  of  the  impact  waves .  Based  on  the  observations  in  our  laboratory,  it  appears  more 
likely  that  large  pressure  and  momentum  fluctuations  in  the  liquid,  that  are  accentuated  near  the 
impingement  point  and  have  their  origins  in  long  wave  instabilities  that  are  present  on  the  jet 
prior  to  impingement,  determine  the  breakup  characteristics  of  turbulent  impinging  jets,  and  are 
responsible  for  the  obvious  differences  between  laminar  and  turbulent  cases.  A  quantitative 
assessment  of  the  specific  mechanisms  controlling  atomization  remains  to  be  established. 

Comparisons  of  the  experimental  results  were  also  made  with  results  from  a  linear 
stability-based  theory.  The  model  predicted  a  monotonically  decreasing  sheet  breakup  length 
with  increasing  Weber  number  which  is  opposite  to  the  trend  observed  for  turbulent  impinging 
jets.  Thus,  the  use  of  the  aerodynamic  instability  model  for  describing  sheet  breakup  appears  to 
be  incorrect,  at  least  for  the  conditions  of  these  studies.  Results  from  computations  from  a  finite- 
difference  Navier-Stokes  model  were  also  compared  to  the  experimental  data  and  interesting 
similarities  between  the  predicted  and  measured  values  of  disturbance  wavelength  were 
observed.  Measurements  of  disturbance  phase  speed  must  be  made  for  further  corroboration  of 
the  model. 

Measurement  of  surface  wave  and  periodic  ligament  formation  phenomena  for  turbulent 
impinging  jets  indicates  that  the  observed  wavelengths  are  directly  proportional  to  the  orifice 
diameter  and  independent  of  jet  velocity  and  impingement  angle,  which  is  also  contrary  to  results 
from  the  linear  stability  model  that  predicts  a  strong  dependence  of  these  variables  on  jet 
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velocity.  It  is  interesting,  however,  that  the  aerodynamic  instability  model  does  provide  quite  a 
reasonable  predictive  capability  for  drop  size.  Measurements  made  in  a  high-pressure  acoustic 
chamber  indicate  that  there  are  some  effects  of  ambient  pressure  at  greater  than  atmospheric 
pressure.  Based  on  the  results  of  the  above  studies,  approaches  to  modeling  impinging  jet 
atomization  should  continue  to  focus  on  pre-impingement  jet  conditions  and  the  physics  near  the 
jet  impingement  point,  specifically  the  source  of  the  observed  impact  waves  that  appear  to 
control  turbulent  impinging  jet  breakup. 

The  experimental  results  were  studied  in  the  context  of  the  Hewitt  Stability  Correlation. 
At  the  conclusion  of  these  cold-flow  studies,  three  observations  with  potentially  high 
significance  can  be  made: 

(1)  The  frequency  with  which  periodic  surface  waves  and  ligament  structures  are  formed 
have  a  marked  similarity  to  the  highest  possible  combustion  instability  frequency  as 
predicted  by  the  stability  correlation. 

(2)  Increases  in  the  value  of  the  stability  parameter,  indicating  increased  stability  margin, 
coincide  with  an  increase  in  measured  mean  drop  size. 

(3)  Increases  in  the  value  of  the  stability  parameter  also  coincide  with  a  increase  in 
polydispersity  of  the  drop  size  distribution. 

Each  of  these  observations  are  consistent  with  a  reasonable  theory  of  combustion  instability. 

Combustion  measurements  at  realistic  chamber  and  injector  conditions  are  critical  to 
defining  combustion  instability  mechanisms.  A  design  of  an  optically-accessible  combustion 
chamber  has  been  developed.  The  first  chamber  will  be  designed  for  250  psi  chamber  pressure. 
Pressure  oscillations  will  be  introduced  into  the  chamber  to  simulate  a  combustion  instability 
event  Injector  elements  will  be  placed  strategically  to  allow  the  discrimination  of  oscillatory 
pressure  effects  from  oscillatory  velocity  effects.  Ethanol  will  be  used  as  the  liquid  fuel. 
Measurements  will  be  phase-locked  to  the  pressure  oscillations,  and  will  include  combustion 
light  measurements  (OH  or  CH  radical  emission)  and  drop  size  and  velocity  measurements.  A 
combustion  response  model  will  be  used  to  analyze  the  combustion  test  results.  Unsteady 
atomization  and  drop  distribution  effects  will  be  accounted  for  in  the  analysis. 
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Abstract 

The  atomization  characteristics  of  sheets  formed  by  both  laminar  and  turbulent  impinging  jets 
were  experimentally  studied  as  a  function  of  flow  and  injector  geometric  parameters.  In  particular,  sheet 
breakup  length  along  the  sheet  centerline,  distance  between  adjacent  waves  apparent  on  the  sheet,  and 
drop  size  distributions  were  measured  over  a  Weber  number  range  between  350  to  6600  and  a  Reynolds 
number  range  between  2800  to  26000.  A  linear  stability-based  model  was  used  to  determine  the  most 
unstable  wavenumber  and  the  corresponding  growth  rate  factor  on  two-dimensional  thinning  inviscid  and 
viscous  sheets.  These  wave  characteristics  were  used  to  predict  both  the  sheet  breakup  length  and  the 
resulting  drop  sizes.  A  second  model,  applicable  for  a  low  Weber  number  regime,  in  which  sheet 
disintegration  is  controlled  by  stationary  antisymmetric  waves,  was  used  to  predict  the  shape  of  the  sheet 
formed  by  two  impinging  liquid  jets.  The  linear  stability-based  theory  predictions  of  breakup  length  did 
not  agree  in  trend  or  magnitude  with  experimental  measurements.  However,  for  Weber  numbers  less 
than  350.  the  measured  breakup  length  for  laminar  impinging  jets  was  within  50%  of  that  predicted  by 
the  stationary  antisymmetric  wave-based  model.  Finally,  drop  size  predictions  based  on  linear  stability 
theory  agreed  in  trend  but  not  in  magnitude  with  the  measured  drop  sizes.  The  contrast  between  the  sheet 
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atomization  characteristics  of  laminar  versus  turbulent  impinging  jets  suggest  that  the  initial  conditions 
of  the  impinging  jets  significantly  influence  the  sheet  breakup  mechanism.  Also,  the  comparison  between 
experimental  results  and  theoretical  predictions  indicates  that  the  impact  wave  generation  process  at  the 
jet  impingement  point  needs  to  be  incorporated  in  the  theoretical  models  for  sheet  atomization. 

Nomenclature 

English  Symbols 
d  diameter 

F  thickness  distribution 

h  sheet  thickness 

k  wavenumber 

/  length 

L  length  of  injection  element 

r  radial  distance  from  impingement  point 

Re  Reynolds  number  ( = U}dJ v,)  based  on  liquid  properties,  jet  velocity,  and  orifice  diamner 

Re,  Reynolds  number  (=U,hlv)  based  on  liquid  properties,  sheet  velocity,  and  sheet  thickness 

s  ratio  of  gas  density  to  liquid  density 

t  time 

U  velocity 

W  maximum  width  of  sheet 

We  Weber  number  (=p,U'dJa)  based  on  liquid  properties,  jet  velocity  and  orifice  Hiamcrcr 

We,  Weber  number  ( =plU,2h/o)  based  on  liquid  properties,  sheet  velocity  and  sheet  thickness 

x  axial  distance  from  impingement  point 

y  coordinate  perpendicular  to  x  in  the  plane  of  the  sheet 


Greek  Symbols 

a  fan  inclination  angle 

0  complex  growth  rate  factor  (0=0r+i0i) 

tj  disturbance  amplitude 

X  wavelength 

n  dynamic  viscosity 

v  kinematic  viscosity 

t  pi  (=3.14159) 

$  angular  coordinate  on  sheet 

p  density 

0  impingement  half-angle 

a  surface  tension 

Subscripts 
10  arithmetic 

30  volume 

b  breakup 

D  drop 

e  edge 

g  gas 

i  imaginary 

j  jet 

/  liquid 

L  ligament 

m  maximum 


A  '  3 


nd  nondimensional 

o  orifice  or  initial 

r  real 

s  sheet 

sw  surface  wave 

Introduction 

Impinging  jet  injectors  are  commonly  used  in  many  rocket  engines,  prominent  examples  being 
the  F-l,  H-I,  Titan  and  XLR-132.1  Furthermore,  laser-drilled  micro-orifice  doublet  impinging  jet 
injectors  have  recently  received  attention  due  to  their  potential  low  cost  and  high  efficiency.1 
Impinging  jet  injectors  are  most  often  used  with  RP-1 /liquid  oxygen  (LOX)  and  nitrogen  tetroxide 
(NTO)/monomethyl  hydrazine  (MMH)  propellant  combinations.3  The  relative  ease  of  fabrication  of 
impinging  jet  injectors  also  makes  this  type  of  injector  an  attractive  alternative  to  coaxial  injector  elements 
typically  used  for  LOX/H,  engines.  Currently,  no  mechanistic  design  analysis  method  exists  for  this 
common  type  of  injector. 

A  schematic  diagram  of  a  typical  like-on-like  impinging  jet  injector  and  the  resultant  spray  is 
shown  in  Fig.  1 .  The  individual  impinging  jet  injector  elements  are  fed  propellant  through  a  manifold. 
The  inlets  are  usually  sharp  edged;  hence  the  flow  detaches  from  the  wall  to  form  a  vena  contraaa. 
The  flow  typically  reattaches  to  the  wall  and  exits  into  the  chamber.  The  two  emerging  cylindrical  jets 
impinge  at  a  point  with  an  impingement  angle  of  20.  which  is  typically  about  60°.  The  Reynolds 
number.  Re  ( =  Udjv,),  and  Weber  number.  We  ( =p,UfdJa ).  based  on  liquid  properties,  jet  velocity,  and 
orifice  diameter  are  on  the  order  of  10*- 10*  under  rocket  conditions.  The  length-to-diameter  ratio  (Udg) 
of  the  orifice  and  the  pre-impingement  length-to-orifice-diameter  ratio  (l/dj  are  both  typically  between 
three  and  ten.  The  symbol  o  shown  in  Fig.  1  is  the  angular  coordinate  in  the  plane  of  the  sheet. 
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The  impinging  liquid  jets  form  a  sheet  in  a  plane  perpendicular  to  the  plane  of  the  jets,  and  the  sheet  is 
usually  canted  at  an  angle,  a,  from  the  normal  to  the  injector  faceplate  to  improve  mixing. 

The  impetus  behind  understanding  the  atomization  mechanism(s)  of  impinging  jet  injectors  is  the 
direct  effect  atomization  has  on  subsequent  combustion  processes,  and  the  established  link  between 
atomization  and  combustion  instability  phenomena.1 3,4  Previous  cold  flow4'*  and  hot-fire3  experimental 
studies  of  impinging  jet  atomization  have  typically  involved  changing  the  injector  geometry  and  flow 
properties  and  recording  the  effect  of  these  changes  on  the  resultant  spray.  Visualizations  of  impinging 
jet  sprays  under  both  cold  flow  and  hot-fire  conditions  indicate  that  the  sheet  breaks  up  in  a  periodic 
manner  into  ligaments,  which  subsequently  disintegrate  into  drops.3  Previous  studies  have  shown  that 
the  sheet  structure,  ligament/drop  formation  frequency  and  drop  size  are  all  sensitive  to  the  injector 
design  and  operating  parameters.4  In  addition,  researchers  have  modeled  the  breakup  behavior  of 
impinging  jets  with  varying  degrees  of  complexity  and  success.6-9  However,  a  fundamental  understanding 
of  the  atomization  mechanising)  still  does  not  exist. 

In  this  study,  measurements  of  sheet  breakup  length,  and  surface  wave  and  periodic  ligament 
separation  characteristics  made  from  instantaneous  images  of  the  impinging  jet  spray,  and  drop  size 
measurements  within  the  spray  field,  are  presented.  These  results  are  compared  to  predictions  made  with 
current  modeling  approaches,  and  with  earlier  experimental  studies.  Similar  to  previous  workers,3 
the  effects  of  laminar  versus  turbulent  jet  conditions  on  the  spray  breakup  process  have  been  studied  over 
a  Weber  number  range  between  350  to  6600  and  a  Reynolds  number  range  between  2800  to  26000. 
The  objectives  of  the  current  study  are  to  extend  and  confirm  previous  experimental  work  on  impinging 
jet  atomization  and  to  provide  a  firm  basis  from  which  to  address  combustion  instability  phenomena 
associated  with  this  type  of  injector.  To  aid  in  the  comparison  of  current  work  with  previous  studies, 
a  review  of  previous  theoretical  and  experimental  work  on  liquid  sheet  breakup  is  presented  next. 


Theoretical  Models  of  Liquid  Sheet  Breakup 


The  study  of  atomization  characteristics  associated  with  many  injector  types  often  involves 
analyzing  the  breakup  behavior  of  liquid  sheets.  There  have  been  a  number  of  theoretical  studies  of  the 
general  case  of  liquid  sheet  atomization. IM0  as  well  as  the  atomization  of  a  sheet  formed  by  two 
impinging  jets.7'9  Models  of  liquid  sheet  atomization  can  be  subdivided  into  numerical 10'12  and  analytical 
treatments13'20  based  on  the  growth  of  infinitesimal  disturbances  due  to  aerodynamic  stresses  on  the  liquid 
sheet  surface.  The  analytical  models  typically  involve  linear  wave  growth,  while  the  numerical  studies 
examine  nonlinear  wave  growth.  In  a  different  approach,  Childs  and  Mansour10  used  a  Navier-Stokes 
method  and  a  Lagrangian  scheme  to  track  the  liquid-gas  interface  to  provide  an  argument  that  boundary 
layer  effects  in  both  the  liquid  and  gas  phases  enhance  the  wave  growth  rate  for  wind  induced 
instabilities.  These  analytical  and  numerical  studies  provide  insight  into  the  important  physical  processes 
causing  liquid  sheet  disintegration,  and  predict  physical  quantities,  such  as  drop  size,  which  compare 
relatively  well  with  observations  made  in  appropriate  experimental  studies.  However,  the  application  of 
linear  and  nonlinear  wave  growth  based  models  to  the  case  of  sheets  formed  by  two  impinging  jets  has 
not  been  thoroughly  investigated.  Anderson  et  al*  adapted  linear  stability  theory  to  the  case  of  impinging 
jet  atomization,  while  Ibrahim  and  Przekwas9  extended  Taylor’s7,19  work  to  predict  the  sheet  shape  at  low 
Weber  numbers.  A  review  of  linear  stability  theory  for  liquid  sheet  atomization  applicable  for  high 
Weber  number  impinging  jets,  as  well  as  Ibrahim  and  Przekwas’9  sheet  shape  model  for  low  Weber 
number  impinging  jets  is  discussed  next. 

Linear  stability  theory  based  on  the  growth  of  infinitesimal  disturbances  due  to  aerodynamic 
stresses  on  the  liquid  sheet  surface  has  been  used  to  describe  the  disintegration  of  liquid  sheets.13'18 
The  disturbance  on  the  sheet  surface.  77,  is  given  by 

-5-  =  ePl'  (1) 
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where  r)0  is  the  initial  displacement  amplitude,  rf,  is  the  growth  rate  and  t  is  time. 


Typically  the  growth  rate.  (3„  is  calculated  for  a  spectrum  of  wavenumbers,  k.  The  disturbance 
wavenumber,  km,  corresponding  to  the  maximum  growth  rate,  fiLm,  controls  the  breakup  process. 
Both  sinuous  (antisymmetric)  and  dilatational  (symmetric)  waves  can  grow;  however,  previous  research 
indicates  that  sinuous  waves  grow  faster  than  dilatational  waves,13  hence  only  the  behavior  of  sinuous 
disturbances  are  considered  hereafter.  The  theory  does  not  predict  a  critical  disturbance  amplitude  for 
sheet  disintegration,  and  consequently,  an  empirical  relation  of  the  following  form  is  used:14 
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where  t„  is  the  sheet  breakup  time,  xb  is  the  breakup  distance  and  Us  is  the  sheet  velocity.  Thus,  through 
the  above  empirical  equation,  the  length  of  the  intact  sheet  can  be  determined  once  a  relation  for  the 
maximum  growth  rate,  is  ascertained.  The  wavenumber,  km,  associated  with  the  fastest  growing 
disturbance  is  used  to  subsequently  predict  the  size  of  ligaments  and  drops  shedding  from  the  edge  of  the 
intact  sheet.14 

The  derivation  of  an  expression  for  the  growth  rate  factor  and  the  wavenumber  has  been  done 
for  an  inviscid  constant  thickness  sheet,13  and  an  attenuating  viscous  sheet.1413  Other  researchers  have 
extended  these  equations  describing  sheet  breakup  to  the  case  of  the  sheets  formed  by  impinging  jets.19 
A  concise  summary  of  the  appropriate  equations  and  their  underlying  assumptions  is  given  next. 

Squire13  investigated  the  growth  of  antisymmetric  disturbances  for  constant  thickness,  inviscid 
liquid  films,  and  derived  expressions  for  the  wavenumber  and  corresponding  growth  rate  coefficient  for 
the  most  unstable  wave.  The  wavenumber,  kn,  for  the  most  unstable  wave  was  given  as 


2a 


where  pt  is  the  gas  density  and  a  is  the  surface  tension.  The  corresponding  maximum  growth  rate  factor 
was  given  by 
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where  p,  is  the  liquid  density  and  h  is  the  sheet  thickness.  The  above  two  relations  are  valid  for 
Wes(=p,Us2h/o)>  1,  and  for  disturbance  wavelengths  large  compared  to  the  sheet  thickness.  Squire13 
compared  his  predictions  of  the  wavelength  of  the  most  unstable  disturbance,  X*  (=2 v/kj,  to  those 
measured  from  photographs  of  liquid  sheets  produced  by  a  nozzle,  and  showed,  in  general,  good 
agreement. 

Dombrowski  and  Johns14  extended  Squire’s13  analysis  by  including  the  effects  of  viscosity  and  a 
diminishing  sheet  thickness  in  analyzing  the  disintegration  of  two-dimensional  liquid  sheets  through 
aerodynamic  instability.  By  considering  a  force  balance  on  the  liquid  between  pressure,  surface  tension, 
viscous  and  inertial  forces,  they  developed  an  equation  relating  the  disturbance  wavenumber  and  the 
growth  rate  factor  to  fluid  and  sheet  parameters. 14  Brodkey21  rewrote  Dombrowski  and  John’s14  general 
dispersion  equation  in  a  clearer  form: 


where 
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where  v,  is  the  liquid  viscosity  and  s  is  the  ratio  of  the  gas  density.to  the  liquid  density.  In  the  above 
equation,  the  first  term  represents  the  rate  of  change  of  momentum  of  a  liquid  element;  the  second  term 
arises  because  of  viscous  effects:  the  third  term  represents  the  effect  of  surface  forces;  and  the  fourth 
term,  the  effect  of  aerodynamic  forces. 


Weihs.13  like  Dombrowski  and  Johns.14  stated  that  Kelvin-Helmholtz  aerodynamic  instabilities 
resulted  in  the  breakup  of  a  thin,  viscous  liquid  sheet.  However,  Weihs13  derived  a  general  solution  for 
the  shape  of  the  sheet  by  employing  hypergeometric  functions,  which  was  subsequently  simplified  to 


apply  to  the  far  field  region  (i.e.  far  from  the  nozzle,  or  the  impingement  point  in  this  case)  where 
changes  in  sheet  thickness  are  negligible.  He  obtained  die  following  equation  for  the  growth  rate  factor: 


The  same  equation  for  the  growth  rate  factor  is  obtained  by  reducing  Dombrowski  and  Johns’14  general 
dispersion  equation  (Eq.  S). 

Equation  6  relates  the  growth  rate  factor,  to  the  disturbance  wavenumber,  k,  for  a  given  fluid 
and  sheet  velocity.  By  differentiating  /3,  with  respect  to  k  and  setting  the  resulting  expression  equal  to 
zero,  the  wavenumber  of  the  fastest  growing  disturbance,  km,  can  be  found.  Once  km  is  known, 
the  corresponding  growth  rate  factor  can  be  calculated  using  Eq.  6. 

To  predict  drop  size  characteristics,  Dombrowski  and  Johns14  reasoned  that  the  sinusoidally- 
shaped  sheet  breaks  into  cylindrical  ligaments  at  crest  and  trough  points  and  related  the  ligament 
diameter,  dL,  to  the  sheet  thickness  and  wavenumber  as: 
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The  ligaments  subsequently  break  up  into  drops  as  a  result  of  surface  tension  induced  symmetrical  wave 
growth.  Dombrowski  and  Johns14  assumed  that  the  wave  grows  until  the  disturbance  amplitude  is  equal 


to  the  ligament  radius,  thus  resulting  in  one  drop  per  wavelength.  The  drop  size,  for  water  in  air  is 
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Taylor7  made  detailed  measurements  of  the  thickness  and  lateral  spread  of  sheets  formed  by 


low-speed  ( <  S.6  m/s)  impinging  water  jets  and  noted7-20  that  the  overall  shape  of  the  sheet  was  related 
to  the  presence  of  stationary  antisymmetric  waves  within  a  limiting  radius,  which  was  defined  as  the 
radius  where  the  Weber  number  based  on  local  sheet  thickness  is  unity.  His  measurements  showed  that 


the  sheet  thickness,  h,  at  any  radial  location,  r.  on  the  sheet  was  of  the  form,  rh=F(d,&),  which  is 
independent  of  Us.  Hasson  and  Peck22  used  Taylor’s7  data  to  verify  their  analytically  derived  expression 
for  the  sheet  thickness. 
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where  d,  is  the  jet  diameter. 

Ibrahim  and  Przekwas9  extended  Taylor’s19  work  to  obtain  an  analytical  solution  for  the  shape 
of  the  sheet  at  low  Weber  numbers  (We<  500),  while  for  the  high  Weber  number  regime  (We >2000), 
they  suggested  the  use  of  Weihs’  analysis.15  In  the  low  Weber  number  regime,  the  authors  suggested  that 
stationary  antisymmetric  waves  determine  the  shape  of  the  sheet,  which  agrees  with  prior  studies  by 
Taylor19  and  Huang.6  In  the  low  Weber  number  regime,  the  expression  for  the  sheet  shape  took  the 
following  form: 
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where  rb  is  the  distance  from  the  impingement  point  to  the  sheet  edge,  and  ha  and  A,  are  the  initial  and 
edge  sheet  thicknesses,  respectively.  Details  for  calculating  the  initial  and  edge  sheet  thicknesses  can  be 
found  in  Ref.  9. 


The  semi-empirical  theoretical  models  discussed  above  predict  the  sheet  shape  for  the  low  Weber 
number  regime,  and  the  sheet  breakup  length  and  drop  size  for  the  high  Weber  number  regime. 
These  parameters  form  the  basis  for  comparisons  with  the  experimental  results  obtained  in  the  present 
study. 


Experimental  Studies  of  Impinging  Liquid  Jets 


Experimental  studies  on  impinging  injector  systems  have  mainly  been  concerned  with  developing 
mean  drop  size  and  mixing  efficiency  correlations  from  cold  flow  measurements.^2*  and  under 
combusting  conditions. 25  Studies  that  shed  light  on  the  important  physical  mechanisms  controlling  sheet 


A  >10 


breakup  and  subsequent  drop  size  formation  have  been  less  numerous;44  important  results  from  these 
studies  are  summarized  below. 

Heidmann  et  al.4  performed  an  extensive  experimental  study  on  the  atomization  characteristics 
of  two  turbulent  impinging  jets  as  a  function  of  orifice  diameter,  de,  length-to-diameter  ratio.  Ud0,  pre¬ 
impingement  length,  ljf  jet  velocity,  U„  impingement  angle,  20,  dynamic  viscosity,  n„  and  surface 
tension,  a.  From  flash  photographs,  Heidmann  et  al.*  identified  four  spray  patterns.  The  first  spray 
pattern  termed  the  closed  rim  regime  was  characterized  by  a  smooth  liquid  sheet  surrounded  by  a  thick 
rim  that  contained  the  major  portion  of  liquid.  This  sheet  pattern  occurred  at  velocities  below  4  m/s. 
The  next  observed  spray  pattern  was  the  periodic  drop  pattern  in  which  waves  were  evident  on  the  sheet 
surface.  In  addition,  drops  detached  tangentially  off  the  sheet  periphery  at  periodic  intervals. 
The  velocity  range  for  this  pattern  was  between  4  to  9  m/s.  The  open  rim  pattern,  also  observed  between 
jet  velocities  of  4  to  9  m/s,  was  characterized  by  a  thinning  sheet,  and  unlike  the  closed  rim  pattern,  the 
outside  rim  did  not  meet  at  the  spray  centerline.  The  last  spray  pattern  identified  was  the  fully  developed 
pattern  in  which  waves  of  drops  were  shed  in  a  periodic  fashion  from  the  sheet  edge.  Fully  developed 
sprays  were  observed  for  jet  velocities  greater  that  10  m/s.  Heidmann  et  al.*  also  noted  that  there  was 
a  sharp  transition  between  open  rim  and  fully  developed  regimes. 

In  addition  to  identifying  the  four  sheet  patterns,  Heidmann  et  al.*  measured  the  shedding 
frequency  of  ligaments  and  drops  from  the  edge  of  the  sheet  and  found  that  this  "wave  frequency”  was 
linearly  proportional  to  (7;  co s0.  The  "wave  frequency"  decreased  with  increasing  impingement  angle, 
while  remaining  relatively  insensitive  to  changes  in  orifice  diameter  and  pre-impingement  length. 
An  important  observation  made  in  this  study  was  the  similarity  between  the  measured  "wave  frequency" 
and  the  frequency  of  high-frequency  instability  modes  observed  in  liquid  rocket  engines.1-4 

Dombrowski  and  Hooper3  performed  an  experimental  study  on  turbulent  and  laminar  impinging 
water  jets.  Dombrowski  and  Hooper3  measured  the  sheet  speed  using  high  speed  cinematography  and 
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drop  size  from  flash  photographs  as  a  function  of  jet  velocity  and  impingement  angle.  Distinct  sheet 
structure  differences  were  seen  between  the  laminar  and  turbulent  impinging  jet  cases.  The  sheet  formed 
from  laminar  impinging  jets  tended  to  produce  much  larger  and  smoother  sheets  as  compared  to  the  sheets 
formed  by  turbulent  impinging  jets.  Measured  sheet  velocities  tended  to  be  between  the  velocity  of  the 
jet  and  the  value  of  Uj  cosd. 

Dombrowski  and  Hooper5  suggested  that  sheet  breakup  does  not  scale  with  Reynolds  number, 
but  is  dependent  on  the  jet  velocity  profile  and  impingement  angle.  For  turbulent  jets,  "impact  waves" 
formed  at  the  impingement  point  disintegrated  the  sheet.  For  the  laminar  impinging  jets,  both  "impact 
waves"  and  aerodynamic  waves  affected  the  sheet  disintegration  process.  The  authors5  also  stated  that 
the  wavelength  of  both  "impact  waves"  and  aerodynamic  waves  and  the  breakup  length  decrease  with 
increasing  jet  velocity.  However,  breakup  length  measurements  made  from  their  photographs  dispute  the 
latter  statement. 

Huang6  experimentally  and  analytically  studied  the  breakup  of  liquid  sheets  formed  by  two 
opposing  water  jets  (28=  180°)  injected  through  standard  ASME  sharp  edged  orifices  at  jet  velocities 
between  2  to  20  m/s.  Huang6  presented  the  measured  nondimensional  breakup  radius  (rJO.SdJ  versus 
Weber  number  ( p,UfdJo ),  and  indicated  two  Weber  number  dependent  breakup  regimes  connected  by 
a  transition  regime. 

In  the  first  regime,  100  <  We<  500,  the  circular  sheets  were  stable  with  a  nearly  perfect  circular 
edge  and  the  nondimensional  ized  breakup  radius  increased  with  increasing  Weber  number. 
Huang6  occasionally  saw  disturbances  originating  from  the  impingement  point  in  this  regime  akin  to  the 
"impact  waves”  observed  by  Dombrowski  and  Hooper.5  Huang6  performed  a  force  balance  between  the 
inertia  and  surface  tension  forces  acting  on  the  circular  sheet  to  obtain  an  expression  that  related  the 
nondimensional  breakup  radius  to  the  Weber  number  (for  We  <900). 
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The  transition  regime  spanned  a  Weber  number  range  from  500  to  2000.  Huang6  noted  that 
cardioid  wa'  ^s  emerged  and  predominated  in  pan  of  this  region  (500  <  We<  800).  The  nondimensionai 
breakup  radius  reached  a  maximum  between  Weber  numbers  of  800  and  1000,  and  antisymmetric  waves 
first  appeared  in  this  regime. 

Huang6  referred  to  the  second  breakup  regime,  We  >  2000  (the  maximum  We  for  his  experiments 
was  30000),  as  the  unstable  liquid  sheet  regime  where  "the  sheet  flaps  with  a  flag-like  motion."  In  this 
regime,  the  nondimensionai  breakup  radius  decreased  with  increasing  Weber  number.  Based  on  an 
analysis  for  a  vibrating  membrane,  a  semi-empirical  relation  for  the  breakup  radius  was  developed6 
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Anderson  et  al .*  investigated  the  spray  characteristics  of  turbulent  impinging  jets  by  measuring 
the  sheet  breakup  length,  x„,  maximum  sheet  width.  W,  and  drop  size  as  a  function  of  flow  velocity  and 
injector  geometry  (26,  d„,  L!d0  and  /,).  The  experimental  apparatus  and  the  operating  conditions  of  the 
impinging  jet  system  were  very  similar  to  those  of  Heidmann  et  al.*  In  this  study,  the  authors*  found  the 
breakup  length  to  increase  with  decreasing  impingement  angle  and  increasing  jet  velocity,  up  to  the 
maximum  velocity  tested  (18.5  m/s),  which  agrees  with  observations  made  in  previous  experimental 
studies.4  The  maximum  sheet  width.  W.  was  relatively  independent  of  jet  velocity  and  impingement 
angle.  Changes  in  orifice  L!d0  did  not  appreciably  affect  the  spray  characteristics,  while  variation  in  the 
pre-impingement  length.  lr  had  a  measurable  effect  on  breakup  length  and  drop  size.  Drop  size  and 
velocity  measurements  made  with  the  Phase  Doppler  Panicle  Analyzer  (PDPA)  at  the  spray  centerline 
downstream  of  the  impingement  point  showed  the  drop  size  decreasing  with  increasing  jet  velocity  and 
increasing  impingement  angle,  and  that  the  drop  velocity  was  nearly  equal  to  the  jet  velocity. 
Analysis  using  linear  stability  theory  provided  predictions  of  drop  size  that  reproduced  the  experimental 
trends:  however,  breakup  length  predictions  did  not  follow  the  experimental  trends. 


Experimental 


The  atomization  characteristics  of  two  impinging  water  jets  were  experimentally  evaluated  using 
the  experimental  arrangement  shown  in  Fig.  2  which  is  similar  to  that  used  by  Heidmann  et  al.4  and 
Dombrowski  and  Hooper.5  and  is  essentially  identical  to  that  used  by  Anderson  et  al .*  The  experimental 
setup  allowed  for  variation  of  the  impingement  angle  and  size  (length  and  diameter)  of  the  precision  bore 
glass  tubes.  Precision  bore  glass  tubes  were  utilized  to  minimize  effects  of  surface  roughness. 
Minor  spatial  adjustments  to  assure  precise  alignment  of  the  impinging  jets  were  provided  by  a 
micrometer  stage  mounted  on  one  of  the  lines.  The  impinging  jets  were  deemed  to  be  aligned  correctly 
when  the  plane  including  the  two  jets  was  visually  observed  to  be  normal  to  the  plane  of  the  resulting 
sheet.  PDPA  measurements  in  the  spray  field  resulted  in  radial  profiles  of  drop  size  that  were  symmetric 
about  the  centerline.8  and  showed  that  this  simple  procedure  yields  correct  alignment.  The  flow  system 
consisted  of  a  9.47-liter  (2.5-gallon)  tank  filled  with  water  pressurized  by  compressed  nitrogen  gas. 
The  flowrate  was  monitored  using  rotameters,  while  the  pressure  drop  across  the  glass  tubes  was 
measured  using  pressure  gauges.  Both  turbulent  and  laminar  impinging  jets  were  studied. 
Experimental  intricacies  of  each  case  are  discussed  below. 

Turbulent  Impinging  Jets 

Precision  bore  glass  tubes  were  attached  to  4.57  mm  inner  diameter,  I.D.,  (6.35  mm  outer 
diameter.  O.D.)  stainless  steel  tubes  with  fittings.  The  diameter  and  length  of  glass  tubes  used  for  the 
turbulent  impinging  jet  experiments  in  this  study  were  the  same  as  those  used  by  Heidmann  et  al* 
The  tube  orifice  diameters  were  0.64  mm.  1.02  mm  and  1.45  mm.  while  the  length  of  the  tubes  was 
50.8  mm:  thus  yielding  a  length-to-diameter  ratio  of  80.  50  and  35.  respectively.  The  inlet  of  the  glass 
tubes  was  sharp  and  the  ratios  of  fitting  diameter  (6.35  mm  l.D.)  to  glass  tube  diameters  were  9.92.  6.23 
and  4.38.  respectively.  The  pre-impingement  length.  /;,  for  these  turbulent  impinging  jet  experiments 
was  set  to  be  25.4  mm.  thus  yielding  pre-impingement  length-to-orifice  diameter  ratios.  ljd0,  of  40,  25 
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and  18  for  the  three  aforementioned  tube  geometries,  respectively.  Measurements  were  made  at 
impingement  angles,  26,  of  40°,  60°  and  80°.  The  velocity  of  the  water  jets  ranged  between  5-20  m/s. 
The  two  water  jets  were  turbulent  since  for  all  the  velocities  and  diameters  studied,  the  Reynolds  number 
was  greater  than  2300.  In  addition,  no  special  contouring  procedures  were  taken  to  ensure  laminar  flow 
and  the  jet  surface  always  appeared  ruffled. 

Instantaneous  images  of  the  spray  field  were  taken  using  diffuse  backlighting  provided  by  a  strobe 
light  with  a  time  duration  of  approximately  5  ps,  while  a  CID  solid  state  camera  (512  by  512  pixels)  was 
used  to  acquire  the  spray  images.  Typically,  17  spray  images  were  obtained  for  each  unique  geometric 
and  operating  condition.  From  the  spray  images,  the  breakup  length,  xb,  and  the  distance  between  any 
adjacent  periodic  structures,  X,  were  measured.  The  breakup  length  was  defined  as  the  distance  from  the 
impingement  point  along  the  spray  centerline  to  where  the  intact  sheet  disintegrates  into  ligaments  and 
drops,  which  is  the  same  definition  given  by  Heidmann  ex  al*  From  the  17  spray  realizations, 
the  average  value  and  corresponding  standard  deviation  of  the  aforementioned  measured  parameters  was 
obtained.  High-velocity  water  jets  emanating  from  small  Uda  tubes  can  lead  to  cavitation.  For  the  cases 
presented  here,  no  cavitation  effects  were  observed.  The  cavitation  number,  defined  as  the  ratio  of  the 
upstream  pressure  less  the  liquid’s  vapor  pressure  to  the  pressure  drop  across  the  orifice,  was  calculated 
for  the  various  measurement  conditions  and  compared  with  Nurick’s  critical  cavitation  number.29 
The  cavitation  behavior  of  the  jets  was  found  to  be  generally  consistent  with  the  results  observed  by 
Nurick29  and  is  more  fully  discussed  in  Ref.  8. 

Drop  size  and  velocity  measurements  were  made  for  the  0.64  mm  diameter  glass  tube  turbulent 
impinging  jet  case  using  an  argon-ion  (X=514.5  nm)  laser-based.  Fast  Fourier  Transform  (FFT)  version 
of  the  Phase  Doppler  Particle  Analyzer  (PDPA).30  31  The  theory  behind  the  operation  of  the  PDPA  has 
been  reviewed  thoroughly  in  the  literature.30-31  hence  only  a  few  system  details  are  mentioned  here. 
The  collection  optics  were  oriented  30°  off  axis  from  the  forward  propagation  direction  of  the  laser  beam. 


which  is  the  optimum  angle  for  measuring  transparent  drops.30  The  optical  configuration  of  the 
transmitting  and  collection  optics  chosen  allowed  the  measurement  drop  size  range  to  span  40  to 
1400  pun. 

Drop  size/velocity  measurements  for  a  set  of  parametric  conditions  were  made  in  a  previous 
study8  at  two  locations  downstream  of  the  impingement  point,  x—  16  and  41  mm.  At  each  downstream 
location,  measurements  were  made  at  6.4  mm  increments  in  the  plane  of  the  sheet  normal  to  the  sheet 
centerline.8  In  this  study,  only  the  drop  size  measurements  at  a  single  spatial  location  (x=  16  mm,  along 
the  sheet  centerline)  are  presented  as  a  function  of  jet  velocity  and  impingement  angle. 
At  each  measurement  location,  8000  drops  were  measured,  since  at  least  5500  data  points  are  needed  for 
±  5%  accuracy  in  mean  diameter  measurements.32 
laminar  Impinging  Jets 

A  series  of  experiments  involving  laminar  impinging  jets  was  also  undertaken  using  the 
experimental  setup  shown  in  Fig.  2.  For  these  tests.  0.51-mm-diameter  precision  bore  glass  tubes  with 
an  UdB  ratio  of  375  were  used.  The  large  Ud0  ratio  ensures  fully  developed  flow.  These  tubes  are  very 
similar  to  those  used  by  Dombrowski  and  Hooper.3  Each  glass  tube  was  attached  to  a  12.7-ram-O.D., 
10.4-mm-I.D.  brass  tube  with  fittings.  To  ensure  laminar  flow  through  the  precision  bore  glass  tubes, 
an  inlet  glass  mbe  was  contoured  such  that  the  fluid  transitioned  smoothly  from  the  brass  tube  to  the  glass 
tube.  The  transition  angle  between  the  inlet  glass  tube  wall  and  the  tube  centerline  was  17°. 
This  contoured  inlet  mbe  was  fused  to  the  0.51-mm-I.D.  precision  bore  glass  mbe. 

Contouring  the  inlet  extended  the  Reynolds  number  range  for  laminar  flow  up  to  10000  as 
indicated  by  pressure  drop  measurements  taken  across  the  glass  tube,  and  the  glass-like  appearance  of  the 
jets.  Often,  a  bursting  phenomenon  ",  attributed  to  velocity  profile  relaxation  effects.33  was  observed 
where  the  smooth  laminar  jet  would  suddenly  and  violently  turn  chaotic  at  some  distance  from  the  mbe 
exit. 
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As  with  the  turbulent  impinging  jets,  instantaneous  images  of  the  laminar  impinging  jet  spray 
were  taken  using  the  CID  camera  and  strobe  light.  Approximately  17  images  at  each  operating  and 
geometric  condition  were  acquired  from  which  the  breakup  length  was  measured.  An  average  value  and 
standard  deviation  of  the  two  measured  quantities  were  determined,  and  compared  to  model  predictions 
and  previous  experimental  studies. 

Results  and  Discussion 

The  disintegration  of  sheets  formed  by  two  impinging  liquid  jets  is  modeled  using  two  existing 
theories;  a  stationary  antisymmetric  wave-based  theory  for  low  Weber  numbers,9  and  a  linear  stability- 
based  theory  for  high  Weber  numbers.8  The  stationary  antisymmetric  wave-based  theory  yielded  the 
shape  of  the  sheet  from  which  the  breakup  length  was  obtained,  while  the  linear  stability-based  theory 
yielded  breakup  length  and  drop  size  predictions.  The  breakup  length  predictions  of  both  theories, 
and  the  drop  size  predictions  from  the  linear  stability-based  model  are  compared  to  experimental 
measurements.  Breakup  length  measurements  were  made  for  both  laminar  and  turbulent  impinging  jets, 
while  drop  size  measurements  were  made  only  for  the  turbulent  impinging  jet  case. 
Finally,  measurements  of  the  distance  between  apparent  wave-like  structures  on  the  sheet  surface  are 
presented  for  the  turbulent  impinging  jet  case  as  a  function  of  jet  velocity  and  orifice  diameter. 

General  Spray  Characteristics 

The  sheets  typically  produced  by  turbulent  impinging  water  jets  are  shown  in  Fig.  3.  Here,  the 
orifice  diameter  was  0.64  mm  and  the  jet  velocities  were  6.4  and  18.5  m/s.  The  free  jet  length  prior  to 
impingement,  or  pre-impingement  length.  lr  was  25  mm,  Ud0  was  80  and  the  impingement  angle,  20, 
was  60°.  The  two  instantaneous  images  shown  in  Fig.  3  fall  into  the  fully  developed  regime  identified 
by  Heidmann  et  al.i  It  is  interesting  to  note  the  apparent  periodic  nature  of  both  the  disturbances  on  the 
sheet  surface  and  the  detaching  ligaments.  Recall  that  Heidmann  et  al .*  measured  the  "wave  frequency" 
of  detaching  ligaments  and  drops  and  observed  it  to  be  linearly  proportional  to  U  cos 0. 


Instantaneous  images  of  sprays  resulting  from  impinging  laminar  water  jets  are  shown  in  Fig.  4. 
In  this  case,  the  orifice  diameter  was  0.51  mm.  the  jet  velocities  were  7.1  and  17.9  m/s,  Ud0  was  375. 
the  impingement  angle  was  60°  and  the  pre-impingement  length,  ljt  was  approximately  10  mm. 
The  images  in  Fig.  4  have  the  same  scale  as  the  images  for  the  turbulent  impinging  jet  case  shown  in 
Fig.  3.  Despite  the  similarity  between  the  operating  and  geometric  parameters  of  the  turbulent  and 
laminar  cases,  the  appearance  of  the  sheets  produced  in  each  case  is  quite  different.  The  sheets  resulting 
from  the  impinging  laminar  jets  are  larger  than  their  turbulent  counterparts.  In  addition,  the  sheets  for 
the  laminar  case  tend  to  be  much  smoother  and  less  chaotic.  Dombrowski  and  Hooper5  also  observed 
distinct  differences  between  turbulent  and  laminar  impinging  jet  cases,  and  attributed  these  differences 
primarily  to  the  different  velocity  profiles  across  the  jets  for  the  two  cases.  Dombrowski  and  Hooper5 
stated  that  "impact  waves"  originating  from  the  impingement  point  were  responsible  for  sheet 
disintegration  for  the  turbulent  impinging  jet  case  as  well  as  for  the  breakup  of  sheets  formed  by  high 
speed  laminar  impinging  jets.  These  "impact  waves"  can  be  seen  in  the  Uj— 17.9  m/s  image  of  Fig.  4, 
and  appear  to  lead  to  the  disintegration  of  the  sheet.  Careful  inspection  of  the  high  velocity  case  in  Fig.  4 
reveals  tears  at  the  center  of  the  intact  sheet  as  well.  In  addition,  the  "impact  waves"  are  dominant  in 
the  center  portion  of  the  sheet.  Dombrowski  and  Hooper5  stated  that  the  slower  moving  fluid  at  the  sheet 
periphery  resulting  from  the  laminar  velocity  profile  across  the  jet  tended  to  damp  out  the  "impact 
waves" . 

The  images  presented  here  demonstrate  the  importance  of  the  jet  condition  prior  to  impingement 
in  understanding  the  atomization  process.  The  differences  between  laminar  and  turbulent  impinging  jet 
conditions  involve  details  related  to  the  velocity  profile  and  turbulence  intensity  of  the  jets  prior  to 
impingement.  Despite  the  clear  importance  of  the  jet’s  initial  condition  on  the  atomization  process, 
no  definitive  studies  presently  exist  to  clearly  differentiate  velocity  profile  and  turbulence  intensity  effects 
on  atomization.  Actual  impinging  injectors  have  small  L!dg  values  (typically  less  than  five)  and  Re 
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numbers  which  lie  in  the  turbulent  regime  (typically  1(f).  Thus,  the  velocity  profiles  developed  under 
these  short  Lld0  ratios  will  not  correspond  to  fully  developed  conditions.  Furthermore,  high  turbulence 
intensity  levels  are  to  be  expected  in  actual  conditions  due  to  the  flow-turning  effects  encountered  in  the 
injector  manifold  of  typical  rocket  engines.  The  degree  to  which  such  phenomena  control  the  atomization 
process  remains  a  serious  challenge  to  experimentalists  and  modelers  interested  in  impinging  jet  injectors. 
Analysis 

The  linear  stability-based  model  was  used  to  predict  breakup  length  and  drop  size.  For  these 
predictions  to  be  made,  the  sheet  velocity,  Us,  must  be  known.  Drop  velocity  measurements  made  in  a 
previous  study  for  turbulent  impinging  jets8  indicated  that  for  the  entire  tested  velocity  range  (S  to  25  m/s) 
and  impingement  angle  range  (40°and  100°),  the  mean  drop  velocity  was  close  to  the  jet  velocity,  Ur 
This  result  suggests  that  the  sheet  speed  is  the  same  as  the  speed  of  the  incoming  jets,  Ur  Taylor7  also 
noted  that  the  jet  and  sheet  speeds  should  be  the  same.  Dombrowski  and  Hooper8  measured  the  sheet 
speed  through  high  speed  cinematography,  and  found  essentially  similar  results  for  the  sheet  velocity. 
Consequently,  the  sheet  velocity,  U„  is  assumed  to  be  the  same  as  the  jet  velocity,  Uj.  In  addition,  it 
is  assumed  that  the  jet  diameter,  djy  is  the  same  as  the  orifice  diameter,  d0. 

The  first  step  taken  was  to  compare  Squire’s13  analysis  with  Dombrowski  and  Johns14  and  Weihsis 
analyses  of  viscous  sheets.  Weihs,13  after  simplifying  his  analysis,  obtained  the  same  expression  relating 
the  growth  rate  factor  and  wavenumber  as  Dombrowski  and  Johns.14  Squire13  provided  explicit 
expressions  for  the  wavenumber,  k„,  (Eq.  3)  and  the  growth  rate  factor,  (3im,  (Eq.  4)  for  the  most 
unstable  wave  which  is  assumed  to  control  the  breakup  process  of  the  liquid  sheet.  For  the  viscous 
sheets,  the  wavenumber  of  the  most  unstable  wave.  was  found  from  Eq.  6,  taken  from  Weihs’ 
analysis.13 

The  wavenumber  and  growth  rate  factor  of  the  fastest  growing  wave  for  both  the  viscous  and 
inviscid  cases  were  calculated  and  compared  as  a  function  of  sheet  velocity  and  sheet  thickness  for  a 
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water  sheet  in  still  air.  With  increasing  sheet  thickness  and  velocity,  the  maximum  wavenumber  predicted 
using  the  inviscid  analysis  becomes  larger  than  that  predicted  by  the  viscous  analysis.  However,  for  the 
cases  studied  here.  Squire’s13  simple  expression  for  the  wavenumber  (Eq.  3)  of  the  most  unstable  wave 
results  in  less  than  1  %  deviation  from  the  viscous  analysis  results. 

Furthermore,  for  the  flow  and  geometric  conditions  of  this  study,  the  maximum  growth  rate 
factor,  {3,  m,  predicted  by  the  inviscid13  and  viscous1415  analyses  were  also  the  same.  However,  for  more 
viscous  liquids  and  higher  gas  densities  (corresponding  to  higher  chamber  pressures),  the  two  analyses 
will  differ.  Based  on  these  observations.  Squire’s  simple  expressions  for  the  wavenumber  (Eq.  3)  and 
growth  rate  factor  (Eq.  4)  of  the  most  unstable  wave  were  used  to  calculate  breakup  length  and  drop  size. 

The  sheet  breakup  length  is  linked  to  the  sheet  velocity,  disturbance  growth  rate,  sheet  thinning 
rate  and  the  ratio  of  the  final  to  the  initial  disturbance  amplitude.  Recall  that  the  breakup  length  is 
defined  as  the  distance  from  the  impingement  point  to  where  the  intact  sheet  disintegrates  along  the  spray 
centerline.  The  empirical  relation  given  by  Eq.  2  is  used  to  compute  the  breakup  length,  jq,,  using  the 
maximum  growth  rate  factor,  0i  m  (see  Eq.  4),  derived  by  Squire.13  Inspection  of  Eq.  4  indicates  the  need 
for  a  sheet  thickness  expression.  For  the  sheet  thickness,  Hasson  and  Peck’s22  expression  (Eq.  9)  was 
used.  With  the  three  aforementioned  equations,  an  explicit  expression  for  the  breakup  length,  xb,  was 
derived  as  follows: 

r  -1  -1 
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where  d0  is  the  orifice  diameter,  s  is  the  ratio  of  gas  density  to  liquid  density.  We  is  the  Weber  number 
based  on  liquid  properties,  jet  velocity  and  orifice  diameter,  and  fi.0)  is  given  as  follows: 

M  -  <13> 

sin3® 

where  6  is  the  impingement  half-angle.  Note  that  the  constant  given  by  Eq.  13  is  contained  in  the  sheet 
thickness  equation  of  Hasson  and  Peck-  (Eq.  9)  for  o=0°  (spray  centerline).  If  a  different  expression 


for  the  sheet  thickness  was  used,  thenfid)  would  have  a  different  form.  For  conditions  involving  more 
viscous  liquids  and/or  higher  gas  densities,  the  expression  for  the  growth  rate  factor,  for  the  viscous 
case  given  by  Eq.  6,  and  Hasson  and  Peck’s22  sheet  thickness  relation  (Eq.  9)  can  also  be  used  in 
conjunction  with  the  empirical  breakup  relation,  Eq.  2.  The  resulting  expression  could  then  be 
numerically  integrated  as  a  function  of  wavenumber  to  solve  for  the  breakup  length.1 
Sheet  Breakup  Length  Measurements 

Sheet  breakup  length  was  calculated  using  Eq.  12,  and  compared  to  breakup  length  measurements 
made  for  both  laminar  and  turbulent  impinging  jets.  The  numerical  integration  procedure  outlined  for 
the  viscous  case  was  also  carried  out  and,  as  expected,  the  breakup  length  was  identical  to  that  predicted 
by  Eq.  12.  The  nondimensional  breakup  length,  xjd0,  is  plotted  as  a  function  of  the  nondimensional 
scaling  parameter  Wl?(l-cos0)2/sin30  in  Fig.  S.  For  the  experimental  measurements,  each  symbol 
represents  an  average  of  17  measurements,  and  the  corresponding  bars  represent  the  ±  standard 
deviations.  The  results  obtained  using  linear  stability-based  theory  collapse  to  a  single  line  when  plotted 
in  the  manner  shown  in  Fig.  5.  However,  the  measurements  do  not  collapse  in  a  similar  fashion. 
Breakup  length  measurements  for  turbulent  impinging  water  jets  (d„=0.64  mm,  Lld0=8O)  were  made  at 
discrete  velocities  of  6.4.  12.4  and  18.5  m/s  and  at  impingement  angles  of  40°,  60°  and  80°.  For  the 
laminar  impinging  jet  case  (d„= 0.51  mm,  Ud0= 375),  breakup  length  measurements  were  made  at 
discrete  velocities  of  7.1,  13.3  and  17.9  m/s,  and  impingement  angles  of  40°,  60°  and  80°.  It  should 
be  noted  that  different  precision  bore  glass  tubes  used  for  the  laminar  cases  sometimes  led  to  different 
breakup  length  results.  This  measurement  variability  was  likely  due  to  small  differences  in  tube  inlets 
(i.e..  transition  angle),  tube  outlets  and  internal  surface  conditions.  The  measurements  were  duplicated 
with  several  sets  of  glass  tubes.  The  sensitivity  of  the  sheet  breakup  process  to  apparently  small 
variations  in  the  tube  characteristics  for  the  laminar  case  reinforce  the  observation  that  the  jet  conditions 
strongly  affect  the  spray  formation  process. 
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The  predicted  nondimensional  breakup  length  decreases  with  an  increase  in  the  scaling  parameter. 
However,  the  opposite  trend  is  observed  for  the  turbulent  case  measurements  for  all  impingement  angles. 
For  the  laminar  impinging  case,  the  nondimensional  breakup  length  increases  to  a  maximum,  then 
decreases  with  increasing  Weber  number.  The  linear  stability-based  model  also  overpredicts  the 
magnitude  of  the  breakup  length.  However,  the  predicted  breakup  length  magnitude  can  be  altered  by 
choosing  a  constant  different  than  the  one  in  Eq.  2. 

Several  conclusions  can  be  drawn  by  comparing  the  laminar  and  turbulent  impinging  jet  cases  to 
each  other  as  well  as  with  the  results  of  previous  studies.  For  the  turbulent  impinging  jet  case, 
the  breakup  length  decreases  with  increasing  impingement  angle,  which  is  opposite  to  the  general  trend 
observed  for  the  laminar  jet  case.  Intuitively,  the  breakup  length  would  be  expected  to  increase  with 
decreasing  impingement  angle  since  an  increasing  amount  of  momentum  is  directed  in  the  axial  direction 
and  the  sheet  is  thicker  on  the  centerline.  Dombrowski  and  Hooper3  indicated  that  ’impact  waves" 
controlled  the  breakup  of  sheets  formed  by  impinging  turbulent  jets,  and  that  the  breakup  length  would 
decrease  with  increasing  jet  velocity;  this  trend  is  not  observed  for  the  cases  studied  here.  The  opposing 
trend  of  breakup  length  as  a  function  of  impingement  angle  for  the  laminar  and  turbulent  cases  suggests 
that  different  breakup  mechanisms  are  operative. 

As  mentioned  previously,  the  breakup  length  increases  to  a  maximum,  then  decreases  for 
increasing  Weber  number  for  the  laminar  impinging  jet  case.  Measurements  showed  that  the  largest 
breakup  length  typically  occurred  between  Weber  numbers  of  SS0-72S.  Huang,6  in  his  study  of  opposed 
impinging  jets  (26=  180°').  showed  that  the  breakup  radius  (or  length)  increased  linearly  to  a  maximum 
value  and  then  decreased  with  increasing  Weber  number.  The  breakup  radius  reached  a  maximum 
between  Weber  numbers  of  800  to  1000.  For  Huang's  case  of  26=  180°,  the  trigonometric  term  in  the 
nondimensional  scaling  parameter.  (l-cos0):/sin30.  is  one:  therefore,  his  data  for  breakup  radius  increases 
linearly  for  the  entire  range  of  the  abscissa  shown  in  Fig.  5.  However,  direct  comparisons  between  his 
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results  and  the  current  data  is  clouded  by  the  different  jet  conditions  of  the  two  sets  of  experiments. 
Huang6  used  sharp-edged  ASME  orifices,  which  indicates  that  the  jet  velocity  profile  was  "plug  flow" 
and  probably  had  low  turbulence  intensity  levels.  Also,  the  equation  derived  by  Huang6  for  breakup 
radius  (Eq.  11)  (We>  1000)  is  similar  to  the  equation  derived  here  (Eq.  12).  Inspection  of  the  two 
equations  shows  that  other  than  the  constants,  the  exponents  of  both  the  We  and  s  terms  are  identical. 
The  constants  differ  because  Huang6  obtained  the  value  for  his  constant  by  curve-fitting  the  equation  to 
the  data. 

Dombrowski  and  Hooper5  suggested  that  "impact  waves"  were  responsible  for  sheet  breakup  for 
intermediate  to  high  speed  laminar  impinging  jets,  and  for  all  but  the  lowest  velocities  for  turbulent 
impinging  jets.  This  idea  seems  plausible  since  the  breakup  length  for  both  the  laminar  and  turbulent 
impinging  jet  cases  approach  similar  values  at  high  jet  velocities  (see  Fig.  5).  In  addition,  flash 
photographs  of  the  sheets  under  the  previously  mentioned  conditions  indicate  the  presence  of  disturbances 
which  appear  to  originate  at  the  impingement  point. 

Ibrahim  and  Przekwas9  proposed  that  their  extension  of  Taylor’s19  stationary  antisymmetric  wave- 
based  theory  should  be  used  in  a  low  Weber  number  regime  (Wie  <500),  while  linear  stability  theory 
should  be  applied  to  a  high  Weber  number  regime  (We>2000).  In  the  low  Weber  number  regime, 
the  shape  of  the  sheet  is  predicted  using  Eq.  10.  The  sheet  shape  predicted  using  the  aforementioned 
equation  is  shown  in  Fig.  6  as  a  function  of  the  impingement  angle  for  an  orifice  diameter  of  0.51  mm 
and  a  jet  velocity  of  7.1  m/s.  Note  that  the  impingement  point  is  at  x=0,  y=0.  The  predicted  sheet 
shape  for  20=60°  looks  similar  to  the  corresponding  image  for  the  same  conditions  shown  in  Fig.  4  (low 
jet  velocity  case.  (7=7.1  m/s).  From  the  shape  of  the  sheet,  the  breakup  length,  x„,  can  be  obtained. 
The  predicted  sheet  shapes  have  pointed  tips,  except  for  the  opposed  jet  case  (20=  180°).  For  the  laminar 
impinging  jet  case  (experimental),  sheets  with  pointed  tips  were  observed  for  all  impingement  angles 


tested  up  to  a  Weber  number  of  about  350.  For  Weber  numbers  greater  than  350,  the  sheets  disintegrated 
before  a  pointed  tip  was  formed. 

Further  perusal  of  Fig.  6  indicates  that  the  breakup  length  increases  with  increasing  impingement 
angle  up  to  about  100°.  after  which  it  decreases.  However,  the  maximum  sheet  width  increases  with 
increasing  impingement  angle  up  to  a  maximum  at  26=  180°.  For  the  impingement  angles  studied  in  the 
present  experiments  (29=40°,  60°  and  80°)  for  the  laminar  impinging  jet  case,  the  general  trend  was  for 
increasing  breakup  length  with  increasing  impingement  angle,  the  same  trend  predicted  by  the  model. 
Since  breakup  length  measurements  were  not  made  for  impingement  angles  greater  than  100°, 
the  decrease  in  breakup  length  with  increasing  impingement  angle  for  impingement  angles  greater  than 
100°  as  predicted  by  the  model  cannot  be  either  confirmed  or  contradicted. 

A  comparison  between  the  breakup  length  predicted  by  linear  stability-based  theory,  Ibrahim  and 
Przekwas’9  extension  of  Taylor’s19  stationary  antisymmetric  wave  theory,  and  the  breakup  length 
measurements  for  the  laminar  and  turbulent  cases  is  shown  in  Fig.  7.  Specifically,  the  nondimensional 
breakup  length,  xjd„,  is  plotted  as  a  function  of  the  non-dimensional  scaling  parameter,  Wefl-costf^/sin3®, 
for  an  impingement  angle  of  60°.  As  seen  in  Fig.  7,  the  breakup  length  as  predicted  by  die  Ibrahim  and 
Przekwas9  model  is  linearly  proportional  to  Weber  number,  whereas  for  the  linear  stability-based  model, 
it  is  proportional  to  We  in.  The  breakup  length  predictions  from  the  Ibrahim  and  Przekwas9  model 
overpredict  the  breakup  length  measurements  made  for  both  the  laminar  and  turbulent  cases;  however, 
the  measured  breakup  length  does  increase  with  increasing  Weber  number,  but  not  as  strongly  as  the 
stationary  antisymmetric  wave-based  model  predicts.  The  model  also  predicts  sheet  shapes  with  pointed 
tips,  as  shown  in  Fig.  6.  These  sheets  were  only  observed  in  the  laminar  impinging  jet  case  for  Weber 
numbers  less  than  350.  which  correspond  to  the  first  two  data  points  for  the  laminar  case  in  Fig.  7. 
As  mentioned  previously,  the  sheets  formed  by  laminar  impinging  jets  disintegrated  before  pointed  tips 
were  formed  for  Weber  numbers  greater  than  350.  thus  explaining  why  longer  sheets  are  predicted  using 
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the  stationary  antisymmetric  wave-based  model.  For  the  turbulent  impinging  jet  case,  the  measured 
nondimensional  breakup  lengths  for  various  jet  diameters  lie  on  nearly  a  single  curve  for  this  fixed 
impingement  angle.  The  same  trend  is  observed  at  the  other  impingement  angles  although  different 
curves  were  observed  for  each  angle.  It  should  be  noted  that  this  collapse  of  the  turbulent  impinging  jet 
results  to  a  single  curve  for  a  fixed  impingement  angle  is  a  reflection  of  the  We  number  dependence  since 
the  geometric  factor  in  the  nondimensional  scaling  factor  remains  constant  for  fixed  impingement  angle. 
Similar  to  Huang’s6  breakup  radius  results  for  opposed  jets,  and  the  present  results  for  laminar  impinging 
jets,  it  is  likely  that  the  breakup  length  for  the  turbulent  impinging  jet  case  will  peak  for  some  Weber 
number  and  then  decrease. 

Drop-Size  Measurements 

Drop-size  measurements  using  the  Phase  Doppler  Particle  Analyzer  (PDPA)  for  the  turbulent 
impinging  jet  cases  were  compared  to  predictions  of  the  linear  stability-based  model.  There  are  two 
models  for  calculating  the  drop  diameter.  The  model  proposed  by  Dombrowsld  and  Johns14  (Eqs.  7  and 
8)  can  be  combined  with  Squire’s13  expression  for  the  fastest  growing  wavenumber  (Eq.  3)  and  the 
nondimensional  breakup  length  expression  derived  earlier  (Eq.  12),  and  recast  in  the  following  form: 

~  =  —x  s'* [We MV  (14) 

(12) 3 

where  dD  is  the  drop  diameter  and  fijS)  is  given  by  Eq.  13.  The  relationship  shows  that  the  ratio  of  the 
drop  diameter  to  the  orifice  diameter  has  a  weak  inverse  dependence  on  the  scaling  parameter,  WeJ[8). 
The  second  model,  proposed  by  Ibrahim  and  Przekwas.4  relates  the  drop  size  to  the  maximum 
wavenumber.  km,  as 

dD  -  f  (15) 

which  can  be  recast  as  follows: 
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(16) 


do  _  2 no  _  2.t 

d°  d0ptUl  s  We 

By  contrasting  the  two  relationships  for  the  nonditnensional  drop  diameter,  it  is  evident  that  the  former 
model  shows  a  dependence  on  the  impingement  angle,  whereas  the  latter  model  is  independent  of  the 
impingement  angle  and  inversely  proportional  to  the  We  number. 

The  drop  size  measurements  are  compared  with  the  linear  stability-based  model  predictions 
(Eq.  14)  in  Fig.  8.  The  drop  size  measurements  shown  in  the  figure  are  for  a  spatial  location,  x=  16  mm, 
along  the  sheet  centerline.  The  abscissa  in  this  figure  is  the  nondimensionai  scaling  parameter.  We  fid), 
whereas  the  ordinate  is  the  drop  diameter  nondimensionai ized  with  the  orifice  diameter.  Note  that  the 
measured  drop  sizes  are  polydispersed.  and  the  arithmetic  mean  diameter,  dl0,  is  taken  for  comparison 
purposes.  The  theory  predicts  a  monodispersed  drop  size  distribution.  The  comparison  shows  that  the 
drop  size  dependence  on  both  the  Weber  number  and  the  impingement  angle  is  similar  for  both 
measurements  and  predictions.  The  nondimensionai izaiion  also  brings  the  data  for  various  impingement 
angles  close  to  a  single  curve.  It  is  not  surprising  that  the  measurements  and  theory  do  not  match 
quantitatively;  better  agreement  would  be  observed  if  a  higher-order  moment  diameter,  say  a  d„,  for  the 
measurements  were  used.  Alternately,  a  larger  empirical  breakup  constant  (currently,  J  /Stwdr=  12) 
would  bring  the  measurements  and  the  model  predictions  closer.  In  this  vein,  an  empirical  breakup 
constant,  J  (3,  „dt =64,  yields  a  semi-empirical  model  of  the  fotm: 

j  •  l*Ls'*[Wefi0)p  (17) 

(64)* 

that  represents  a  least-squares  fit  to  the  measured  dl0,  as  shown  in  Fig.  8.  Finally,  the  use  of  the  second 
drop-size  model  (Eq.  16)  predicts  drop  sizes  which  are  an  order  of  magnitude  greater  than  those 
measured.  Therefore,  based  on  these  observations,  the  former  model  for  drop  size  has  more  merit. 
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Surface  Wave  Measurements 


Heidmann  et  al .*  measured  the  "wave  frequency"  of  detaching  ligaments  and  drops  from  the  edge 
of  the  intact  sheets  formed  by  impinging  jets  as  a  function  of  operating  and  injector  parameters. 
Two  important  results  of  their  study4  were  that  the  "wave  frequency"  was  linearly  proportional  to  U1  cos0. 
and  the  measured  "wave  frequency  "  closely  matched  the  frequencies  associated  with  observed  combustion 
instability  phenomena. 

Because  of  the  interesting  results  and  the  potential  implications  of  the  Heidmann  et  al.*  study, 
similar  measurements  were  made  in  this  study.  Specifically,  for  the  turbulent  impinging  jet  case, 
the  distance  between  apparent  adjacent  periodic  structures  on  and  beyond  the  intact  sheet  surface  were 
measured  as  a  function  of  flow  and  injector  characteristics.  The  types  of  periodic  structures  measured 
include  surface  waves,  edge  ligaments  and  detached  ligaments,  all  of  which  are  shown  in  Fig.  9. 

The  surface  wave  category  included  disturbances  on  the  intact  sheet  surface.  Typically,  the 
distance  between  adjacent  surface  waves,  increased  with  increasing  distance  from  the  impingement 
point.  Separation  distance  measurements  were  made  over  the  entile  intact  sheet  whenever  possible  to 
obtain  an  average  Xw.  The  separation  distance  of  adjacent  surface  waves  is  plotted  as  a  function  of  jet 
velocity  and  orifice  diameter  for  an  impingement  angle  of  60°  in  Fig.  10.  Each  point  plotted  represents 
an  average  between  10  and  50  separate  measurements,  depending  on  operating  conditions,  while  the 
corresponding  bars  represent  the  ±  standard  deviation.  Despite  the  large  standard  deviation,  there  is  a 
distinct  increase  in  with  increasing  orifice  diameter.  However,  for  the  cases  studied,  the  separation 
distance  was  relatively  insensitive  to  changes  in  jet  velocity.  Although  not  presented  in  graphical  format, 
it  was  also  found  that  the  separation  distance,  Xw,  was  insensitive  to  changes  in  impingement  angle. 

Similar  measurements  were  made  for  the  edge  and  detached  ligament  categories.  Edge  ligaments 
were  defined  as  strands  of  liquid  attached  to  the  intact  sheet  periphery.  Detached  ligaments  were  strands 
of  fluid  either  completely  or  nearly  completely  separated  from  the  downstream  edge  of  the  intact  sheet. 
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In  general,  the  average  separation  distance  between  adjacent  detached  ligaments  was  larger  than  that  of 
the  edge  ligaments.  Likewise,  the  average  separation  distance  between  adjacent  edge  ligaments  was  larger 
than  that  of  the  surface  waves.  Again,  the  trends  observed  for  were  also  observed  for  the  edge 
ligament  category.  The  behavior  of  the  separation  distance  between  adjacent  detached  ligaments  was  not 
as  apparent  as  the  other  two  categories  since  not  enough  measurements  were  made  to  yield  a  meaningful 
average  and  standard  deviation.  It  is  interesting  to  note  that  the  measured  separation  distances  of  the 
surface  waves  and  edge  ligaments  are  similar  in  magnitude  to  the  fastest  growing  wavelength  predicted 
by  Squire13 

4.  =  —2  (18) 

rtf 

However,  the  wavelength  of  the  most  unstable  wave  predicted  by  Eq.  18  is  strongly  dependent  on  sheet 
velocity  and  independent  of  the  orifice  diameter.  This  is  contrary  to  the  observed  behavior  of  the  various 
measured  separation  distances. 

Summary  and  Conclusions 

A  systematic  study  of  the  atomization  of  impinging  liquid  jets  investigating  the  effects  of  jet 
conditions  (laminar  versus  turbulent),  orifice  diameter,  impingement  angle  and  jet  velocity  has  been 
conducted.  Results  of  the  present  study  have  been  compared  to  current  theories  in  terms  of  sheet  breakup 
length,  drop  size  and  sheet  shape.  Experiments  contrasting  laminar  and  turbulent  jet  conditions  clearly 
demonstrate  that  the  jet  conditions  have  a  dramatic  effect  on  the  atomization  process.  Specifically,  the 
measured  breakup  lengths  for  the  laminar  impinging  jet  case  are  longer  and  displayed  different  trends  as 
a  function  of  jet  velocity  and  impingement  angle  than  those  of  the  turbulent  impinging  jet  case.  The 
present  results  are  in  complete  agreement  with  earlier  studies  by  Dombrowski  and  Hooper3  with  regards 
to  the  effects  of  jet  conditions  on  impinging  jet  atomization.  These  results  encourage  speculation  that 
velocity  profile  and  turbulence  characteristics  of  the  jets  strongly  affect  the  atomization  processes  for 
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impinging  jet  injectors.  However,  quantitative  assessment  of  the  specific  mechanisms  controlling 
atomization  remain  to  be  established. 

For  low  Weber  number  (We <350)  laminar  impinging  jet  conditions,  predictions  based  on 
stationary  antisymmetric  wave-based  theory  were  within  50%  of  the  observed  breakup  length. 
Comparisons  were  also  made  with  a  linear  stability-based  theory.  The  linear  stability-based  theory 
predicted  a  monotonically  decreasing  sheet  breakup  length  with  increasing  Weber  number  which  is 
opposite  to  the  trend  observed  for  turbulent  impinging  jets.  In  fact,  experimental  observations  for  both 
the  laminar  and  turbulent  jets  argue  strongly  for  the  "impact  wave"  theory  put  forward  by  Dombrowski 
and  Hooper3  as  the  operative  mechanism  leading  to  sheet  breakup.  Thus,  the  current  use  of  linear 
stability-based  theories  for  describing  sheet  breakup  for  impinging  jet  conditions  appears  to  be 
questionable  for  the  range  of  conditions  investigated  in  this  study.  In  contrast,  the  linear  stability-based 
theory  provided  reasonable  predictive  capability  for  drop  size  with  respect  to  the  trends  observed  for 
increasing  Weber  number  and  impingement  angle.  Measurements  regarding  the  surface  wave  and 
periodic  ligament  formation  for  the  turbulent  jet  studies  indicate  that  the  observed  wavelengths  are  directly 
proportional  to  the  orifice  diameter  and  independent  of  jet  velocity  and  impingement  angle. 

Based  on  the  results  of  the  present  study,  approaches  to  modeling  impinging  jet  atomization 
should  focus  on  including  the  "impact  wave"  process  identified  by  Dombrowski  and  Hooper3,  which  none 
of  the  current  models  treat.  Periodic  perturbations  associated  with  the  jet  inertia  at  the  impingement 
region  are  often  argued  to  be  the  likely  mechanism  for  generating  "impact  waves".  A  fundamental 
mechanistic  model  for  the  generation  and  growth  of  "impact  waves"  and  their  association  with  the 
subsequent  atomization  process  is  currently  lacking.  Additionally,  the  relationship  between  the  periodic 
surface  wave  and  ligament  structures  observed  in  the  present  work  needs  to  be  considered  in  light  of  their 
potential  to  initiate  and  sustain  combustion  instability  phenomena  associated  with  impinging  jet  rocket 
injectors. 
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Figure  Captions 

Fig.  1.  Schematic  diagram  of  a  typical  impinging  jet  spray. 

Fig.  2.  Schematic  diagram  of  experimental  arrangement,  a)  top  view;  b)  side  view. 

Fig.  3.  Instantaneous  images  of  sprays  formed  by  two  turbulem  impinging  water  jets.  The  jets  issued 
from  0.64  mm  inner  diameter,  Lld0- 80  precision  bore  glass  tubes.  The  impingement  angle,  28,  was 
60°,  while  the  pre-impingement  length.  ljt  was  25  mm.  Flow  is  from  left  to  right. 

Fig.  4.  Instantaneous  images  of  sprays  formed  by  two  laminar  impinging  water  jets.  The  jets  emanated 
from  0.51  mm  inner  diameter,  L/d0=375  precision  bore  glass  tubes.  The  impingement  angle,  28,  was 
60° ,  while  the  pre-impingement  length,  lr  was  approximately  10  mm.  Note  that  the  scale  of  these  images 
is  the  same  as  those  in  Fig.  3.  In  addition,  the  orifice  diameter  and  jet  velocities  are  similar  for  the 
images  in  both  Fig.  3  and  this  figure.  Flow  is  from  left  to  right. 

Fig.  5.  Nondimensional  breakup  length,  xjd0,  plotted  as  a  function  of  the  scaling  parameter 
We(l-cos0)2/sin30,  where  We  is  the  Weber  number  and  8  is  the  impingement  half-angle.  The  solid  line 
represents  predictions  made  using  linear  stability-based  theory.  The  solid  symbols  represent  the  turbulent 
impinging  jet  case  measurements  (da= 0.64  mm),  while  the  hollow  symbols  represent  the  laminar 
impinging  jet  case  measurements  (d„=0.51  mm).  Each  point  plotted  represents  the  mean  value  of  17 
separate  measurements,  and  the  bars  indicate  the  ±  standard  deviation  of  the  measurements  for  each 
measurement  condition. 

Fig.  6.  Sheet  shape  predictions  for  water  in  air  from  the  Ibrahim  and  Przekwas9  model  plotted  as  a 
function  of  impingement  angle.  The  orifice  diameter  was  0.51  mm.  and  the  jet  velocity  was  7.1  m/s. 
The  predicted  sheet  shape  for  the  60°  impingement  angle  case  should  be  compared  to  the  corresponding 
spray  image  (Fig.  4a).  Note  that  the  scales  are  different. 
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Fig.  7.  Nondimensional  breakup  length,  xjda,  plotted  as  a  function  of  the  scaling  parameter 
WWl-costf^/sin3#.  where  We  is  the  Weber  number  and  6  is  the  impingement  half-angle.  The  solid  line 
represents  predictions  made  using  linear  stability-based  theory.  The  broken  line  represents  breakup  length 
predictions  made  using  the  stationary  antisymmetric  wave-based  model  for  an  impingement  angle  of  60°. 
The  solid  symbols  represent  the  turbulent  impinging  jet  case  measurements,  while  the  hollow  symbols 
represent  the  laminar  impinging  jet  case  measurements.  The  experimental  measurements  shown  are  for 
an  impingement  angle,  26,  of  60°.  Each  point  plotted  represents  the  mean  value  of  17  separate 
measurements,  and  the  bars  indicate  the  ±  standard  deviation  of  the  measurements  for  each  measurement 
condition. 

Fig.  8.  Nondimensional  drop  size,  djd0,  plotted  as  a  function  of  a  nondimensional  scaling  factor  and 
impingement  angle.  Drop  size  predictions  made  using  linear  stability-based  theory  and  drop  size 
measurements  made  using  the  PDPA  at  an  axial  location  of  x=  16  mm,  along  the  spray  centerline  (y=0), 
are  shown. 

Fig.  9.  A  typical  spray  formed  by  two  turbulent  impinging  water  jets  issuing  from  0.64  mm  diameter, 
Lld0= 80  glass  tubes.  The  impingement  angle  was  60°  and  the  pre-impingement  length  was  25  mm. 
Flow  is  from  left  to  right.  The  distance  between  adjacent  surface  waves,  edge  ligaments  and  detached 
ligaments  were  measured  from  such  spray  images. 

Fig.  10.  The  measured  separation  distance  between  adjacent  surface  waves,  plotted  as  a  function  of 
jet  velocity  and  orifice  diameter  for  the  turbulent  impinging  jet  case  at  an  impingement  angle  of  60° . 
Each  point  plotted  represents  an  average  between  10  to  50  separate  measurements,  while  the 
corresponding  bars  indicate  the  associated  standard  deviation. 
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